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ABSTRACT 


The  Fifth  Technical  Report  is  a  continuation  of  and  an  addition 
to  the  First  (Progress),  Second,  Third,  and  Fourth  Technical  Reports 
to  the  Office  of  Naval  Research,  dated  April  30,  1968,  July  15,  1968,  and 
August  31,  1971  respectively. 

The  realization  that  in  mechanical  components  stress,  strength, 
and  the  factors  that  modify  them  arc  by  nature  distributed  has  pointed 
to  the  need  for  design  methods  that  incorporate  the  unique  powers  of 
the  thoeries  of  probability  and  statistics.  Furthermore,  since  with 
any  designed  product  of  a  mechanical  nature  a  finite  probability  of 
failure  is  associated,  this  suggests  that  another  probability.  Reli¬ 
ability,  must  be  the  logical  measure  of  design  adequacy  of  such 
products . 

As  the  design  of  reliable  and  optimized  components  and  systems 
of  dynamic  and  rotary  machinery  is  becoming  a  matter  of  increasing 
concern,  methods  by  which  a  specified  level  of  reliability  can,  in  an 
priori  manner,  be  designed  into  a  component  assume  great  importance. 

So  also  does  the  existence  of  supporting  design  data  in  distributional 
form. 

The  Fifth  Technical  Report  describes  the  continued  experimental 
research  on  material  behavior  under  fatigue  testing.  To  assure  the 
accuracy  and  validity  of  results,  the  calibration  of  the  two  Wiedemann 
Reversed  Bending  Fatigue  Reliability  Research  Machines  and  the  four 
Wire  Fatigue  Reliability  Research  Machines  was  reaccomplished. 
Descriptions  and  results  of  the  calibration  are  presented.  Cycles-to- 


nz 


failure  testing  was  accomplished  on  fine  wire,  0.040  inch  diameter, 

made  of  AISI  4130,  1058,  1018  and  4540  steel.  The- Wiedemann. Machines 

were  used  to  perforin  cycles-to-failure  testing  of  AISI  4130  steel 

grooved  rod,  0.0937  inch  base  diameter.  Also,  endurance  strength, 

staircase  tests  were  accomplished  on  AISI  1038  steel  grooved  rods  of 

0.2700  inch  base  diameter.  These  endurance  strength  results  were  used 

in  a  notch  sensitivity  analysis  described  herein.  Axial  fatigue 

testing  was  done  on  AISI  1018  steel  grooved  rods  to  generate  a  Goodman 

strength  diagram.  Comparison  of  the  experimental  diagram  with  two 

theoretical  diagrams  is  described.  The  results  of  all  testing  were 

analyzed  to  determine  the  distributional  parameters.  The  analysis 

techniques  which  fit  the  normal,  log  normal  and  Weibull  distributions 

c 

to  the  data  are  summarized  and  the  computer  programs  are  included  in 
the  Appendices.  Individual  test  results  are  listed  in  the  data  section 
and  the  results  of  the  analysis  of  distributions  are  shown  in  histograms 
overscribed  with  the  frequency  distribution  curves.  Interpretation  of 
the  tests  results,  conclusions  drawn  therefrom,  and  recommendations 
which  identify  areas  where  further  research  is  needed,  are  included. 
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iml'tin  * : 


1.0  INTRODUCTION 


1.1  PRIOR  RESEARCH 


Research  under  the  Office  of  Naval  Research  Contract  N00014-67- 
A-0209-002  was  initiated  on  1  February  1967.  This  program  has  in¬ 
cluded  both  theoretical  and  experimental  research.  In  the  -theoretical 
research,  the  interactions  among  stress  and  strength  phenomena  in  the 
design  and  life  of  dynamic  and  rotary  machinery  have  been  investigated, 
and  a  method  was  developed  for  applying  the  strength-stress  relation¬ 
ship  to  optimize  the  design  of  machinery  for  a  finite  life  at  a 
specified  level  of  reliability.  The  experimental  research  program  has 
included  the  testing  of  specimens  of  various  materials  and  geometry  to 
support  and  validate  the  results  of  the  theoretical  research. 

During  the  first  reporting  period  a  literature  search  was 
conducted  to  survey  the  extent  and  nature  of  the  interactions  among 
the  design  phenomena  important  to  dynamic  and  rotary  machinery,  to 
assess  the  availability  of  data  in  statistical  form,  and  to  assess  the 
state  of  development  of  the  probabilistic  design-by-reliability  theory 
and  methods.  The  Design-by-Reliability  Methodology  was  discussed  in 
the  First  Progress  Report  (1),  and  the  applicable  mathematical 
equations  were  given.  That  report  also  included  discussions  on  the 
definition  of  failure  mechanisms  in  shafts,  bearings,  and  other 
mechanical  components  and  the  methodology  whereby  conventional  designs 
could  be  improved.  An  experimental  research  program  was  proposed  for 
the  development  of  statistically  described  fatigue  strength  properties 
of  materials.  Testing  machines  to  accomplish  the  experimental  program 
were  identified  and  their  availability  was  confirmed. 

During  the  second  reporting  period,  1  February  1968  -  31  January 
1969,  the  testing  machines  were  calibrated  and  experimental  fatigue 
testing  was  initiated.  The  initial  material  selected  for  study  was 
AJSI  4340  cold  rolled  and  annealed  steel.  Specimens  appropriate  for 
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each  type  of  testing  machine  were  obtained  and  their  physical  and 
geometric  properties  were  evaluated.  The  description  of  the  machine 
calibration  procedures  and  results,  the  procurement  and  preparation  of 
test  specimens  and  the  results  of  fatigue  testing  accomplished  during 
the  second  period  are  described  in  the  Second  Technical  Report  (2). 

During  the  third  reporting  period,  1  February  1969  -  31  August 
1970,  emphasis  was  placed  upon  generating  statistical  data  needed  in 
the  design  of  dynamic  and  rotary  machinery  for  specified  levels  of 
reliability.  Static  tensile,  and  endurance  strength  parameter  dis¬ 
tributions  and  cycles-to-failure  distributions  at  fixed  stress  levels 
were  obtained  for  AISI  4340  and  4130  steel.  Specimens  tested  included 
0.0625  in.  diameter  wires  subjected  to  reversed  bending  on  the  Wire 
Fatigue  Machines,  0.0937  in.  diameter  rods  subjected  to  reversed 
bending  on  the  Wiedemann  Machine,  0.065  in.  diameter  rods  subjected  to 
combined  mean  and  alternating  axial  loads  on  the  Axial  Fatigue  Machine, 
and  0.270  in.  diameter  rods  subjected  to  reversed  bending  on  the  Ann 
Arbor  Machine.  The  Wire  Fatigue  Machines  were  recalibrated  for  the 
diameter  of  wires  being  tested,  and  the  Ann  Arbor  Machine  was  cali¬ 
brated.  A  literature  search  was  conducted  to  compile  the  published 
data  on  the  materials  tested  and  the  results  of  the  experimental 
program  were  compared  with  the  published  data  as  part  of  the  analysis 
program.  Distributional  S-N  diagrams  were  prepared  which  are  directly 
useable  by  designers.  Descriptions  and  results  of  the  research 
accomplished  during  the  third  operating  period  is  contained  in  the 
Third  Technical  Report  (3) . 

During  the  period  1  February  1970  -  31  August  1971  the  exper¬ 
imental  research  was  expanded  to  include  specimens  made  of  two  new 
materials:  AISI  1018  and  AISI  1038  steel.  Testing  included  static 
tensile  tests,  and  hardness  and  geometry  measurements  needed  to 
determine  the  distributional  parameters  of  strength  and  dimensions. 
Fatigue  tests  included  both  endurance-strength  tests  and  cycles-to- 
failure  tests  needed  for  the  preparation  of  Goodman  strength  diagrams 
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and  S-N  diagrams.  Concurrent  research  included  an  investigation  of 
stress  concentration  through  photoelastic  models.  The  analysis 
verified  that  specimens  fabricated  with  relief  angles  and  subjected 
to  fatigue  tests  provided  valid  data.  In  order  to  strengthen  the 
analysis  of  fatigue  testing  data,  a  manual  method  of  fitting  a  three- 
parameter  Weibull  distribution  was  developed.  Thus  the  evaluation  of 
the  distributions  of  fatigue  data  would  include  the  normal,  the  log¬ 
normal,  and  the  Weibull  distributions.  In  addition  to  the  usual 
goodness -of-f it  tests  for  the  normal  and  lognormal  distributions,  an 
automatic  plot  of  a  data  histogram  overscribed  by  a  theoretical  curve 
having  the  estimated  parameters  was  developed  and  used.  Specific 
examples  using  the  data  generated  with  design-by-reliability  method 
to  optimize  designs  of  rotating  mechanical  components  were  developed 
and  included  in  the  Fourth  Technical  Report  (4) . 

1.2  OBJECTIVES  FOR  THIS  REPORT 

The  objectives  for  the  period  1  September  1971  -  31  August  1972 
were  to  continue  the  theoretical  and  experimental  research  of  wire  and 
rod  steel  specimens  in  order  to  provide  designers  with  the  capability 
to  improve  the  design  of  dynamic  rotary  mechanical  components.  Specific 
objectives  for  each  type  of  research  machine  were  as  follows: 

1.2.1  WIRE  RESEARCH  MACHINES 

The  objective  was  to  generate  static  and  dynamic 
strength  distributions  for  wire  specimens  and  to  compile  the  data  in 
a  form  directly  useable  by  designers.  Tasks  conducted  in  order  to 
achieve  the  objective  included: 

1.  Static  tensile  tests  to  determine  the  distributions  of  the 
yield, ultimate  and  breaking  strengths  of  wire  specimens  of 
AISI  4340,  4130,  1018  and  1038  steels. 


4 


2.  Cycles -to -failure  tests  to  determine  the  distribution 
of  fatigue  cycle  life  at  various  levels  of  reversed 
alternating  bending  stress. 

3.  Stress-to-failure  tests  to  determine  the  distribution 
of  fatigue  endurance  strength  for  3  x  10^  cycles. 

4.  Construction  of  distributional  S-N  diagrams  for  the 
given  specimen  materials  and  geometry. 


1.2.2  WIEDEMANN  RESEARCH  MACHINES 

The  objective  was  to  generate  and  analyze  dynamic 
distributional  strength  data  for  AISI  4130  and  1058  steel  rods  with 
various  groove  radii. 

Tasks  included: 

1.  Cycles-to-failure  tests  to  determine  the  distribution  of 
fatigue  life  for  each  specimen  material  and  geometry  sub¬ 
jected  to  reversed  alternating  bending  stress. 

2.  Stress-to-failure  tests  to  determine  the  distribution  of 
the  endurance  strength  for  3  x  10^  cycles  of  life. 

3.  Construction  of  distributional  S-N  diagrams  for  the  specimen 
materials  and  geometries  tested. 

4.  Analyze  the  endurance  strength  parameters  for  the  different 
specimen  groove  radii  to  determine  the  effect  of  groove 
radii  on  fatigue  life. 

5.  Develop  notch  sensitivity  and  fatigue  stress  concentration 
factors  for  the  different  specimen  materials  and  geometries 
tested. 

6.  Completion  of  a  theoretical  study  of  notch  sensitivity  and 
associated  stress  concentration  factors  and  comparison  with 
the  distributional  values  obtained  from  the  experimental 
data. 


a.  *■***- * 


1.2.3  AXIAL  FATIGUE  RESEARCH  MACHINE 


The  objective  via"  j  generate  endurance  strength  dis¬ 
tributions  for  AISI  1018  steel  rods  and  develop  the  associated 
distributional  Goodman  diagram. 

Tasks  included: 

1.  Staircase  tests  to  determine  the  distribution  of  the 
endurance  strength  for  2  x  10  cycles  of  life  for 
specimens  subjected  to  various  ratios  of  axial  alter¬ 
nating  and  mean  loads. 

2.  Construction  of  a  distributional  Goodman  strength 
diagram. 

3.  Determination  of  failure  governing  criteria  which 
correlate  with  the  Goodman  strength  diagram. 
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2.0  EXPERIMENTAL  RESEARCH  PROGRAM  AND  DATA 

The  experimental  program,  undertaken  during  the  Fifth  Phase  of 
this  research  program  is  presented  in  Sections  2.1  through  2.3.  A 
cross-reference  between  the  test  program  and  the  generated  and  reduced 
data  is  also  presented.  The  static,  cycles- to-failure,  and  staircase 
strength  data  was  used  to  generate  S-N  diagrams,  Goodman  diagrams,  and 
fatigue  stress  concentration  and  notch  sensitivity  factors. 

Table  2.0-1  provides  a  list  of  the  research  completed,  in 
progress,  or  scheduled  for  initiation.  In  the  next  final  report  a 
cross-reference  of  all  of  the  research  completed  will  be  provided. 

This  table  will  be  broken  down  by  Data  State,  i.e..  Raw  Data,  Reduced 
Data,  Staircase  Plot,  and  Distribution  Plot,  and  by  Data  Type,  i.e.. 
Static  Ultimate  Strength,  Staircase,  Endurance,  Break  Length,  Diameter, 
Hardness,  Cycles  to  Failure,  Yield  Strength,  and  Elongation. 

2.1  WIRE  RESEARCH  MACHINE 

The  wire  research  machines  were  the  same  as  those  used  in 
previous  phases  of  this  research.  The  calibration,  and  the  details 
about  the  research  program  and  data  for  this  period  of  research  are 
given  next. 

2.1.1  CALIBRATION  OF  THE  WIRE  FATIGUE  RESEARCH  MACHINES 

The  calibration  procedures  for  the  wire  fatigue 
rest  arch  machines  were  the  same  as  those  described  in  the  Third 
Tet.'.i«,:cal  Report  [3,  pp.  38-84],  The  calibration  data  for  Machines 
Number  1,  2,  3,  and  4  appear  in  Tables  10.1-1  through  10.1-4  of 
Section  II  of  this  report. 


Preceding  page  blank 


Test  complete 
Test  in  progress 
Test  planned 


Table  2.0-1  <Cont„)  Summary  of  experimental  research- program  fo-r  the  Office  of  Naval  Research 
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Table  2.0-1  (Cent.)  Summary  of  experimental  research  program  for  the  Office  of  Naval  Research 
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The  static  stress  versus  measured  strain  data  are  given  in 
Table  8- 1-1)  and  summarized  in  Table  9.1-2.  Hie  calculations  for 
the  standard  deviation  of  strain  are  given  in  Table  9.1-3.  The 
results  of  the  static  axial  calculations  which  related  the  measured 
strain  to  a  known  stress  are  presented  in  Table  9.1-4.  A  "Stress 
Versus  Measured  Strain"  chart  was  constructed  from  the  values  in 
Table  9.1-4  and  is  given  in  Figure  9.1-1  .  The  line  in  Figure 
9  .1-1  gives  the  relationship  between  the  mean  strain  and  pan  weight. 
This  also  verifies  the  linearity  of  the  strain  gages  used.  The  least 
squares  technique  was  used  to  calculate  the  best  fit  straight  line 
giving  a  Y  (stress)  intercept  of  4,267  psi  and  a  slope  of  55,379  x  106 
psi/inch/inch  of  strain. 

Thus  the  equation  which  relates  the  measured  strain  in  p-in./in. 
to  the  actual  stress  in  psi  in  the  0.040  in.  diameter  wire  specimens 
is: 

Stress  =  Measured  strain  x  55,379  x  106  +  4,267.  (2.1-1) 

Eq.  (2.1. -1)  provides  the  means  for  determining  the  applied 
stress  for  each  deflection  angle  utilized  in  the  in-machine  calibration. 
The  mean  strain  measured  for  the  selected  deflection  angles  for  all 
machines,  shown  in  Table  9.1-1  ,  were  converted  to  stresses  as  shown 
in  Table  9. 1-5  and  plotted  for  each  machine  in  Figs.  2.1-1  through 
2.1-4  to  provide  the  calibration  curves.  The  curves  compared  well 
with  the  results  of  the  previous  calibration. 

2.1.2  RESEARCH  PROGRAM  AND  DATA 

The  Table  2.1-1  provides  a  cross-reference  between  the 
test  program  and  the  data  generated  and  reduced.  The  raw  and  reduced 
data  can  be  found  in  Volume  II  of  this  report.  Each  data  type  is 
broken  down  by  Wire  Fatigue  Machines,  Wiedemann  Fatigue  Machines,  and 
Axial  Fatigue  Machine,  and  by  Materials  AISI  4130,  4340,  1038  and 
1018  steel. 
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Carriage  Setting 


Using  the  Circular  Arc  Scale  in  Degrees 


Fig*  2.1-1  Calibration  chart  for  Wire 


Machine  #1 . 
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0.040"  Diarceter  Wire-AISI  1038  Steel. 


Carriage  Setting  Using  the. Circular  Arc  Scale  in  Degre* 


Flg*  2.1-2  Calibration  chart  for  Wire  Machine  #2. 


Maximum  Bending  Tensile  Stress  in  kpsi 


Test  Program  for  Wire  Fatigue  Machines 
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2.2  WIEDEMANN  RESEARCH  MACHINES 

The  research  machines  were  the  same  as  those  used  in  previous 
phases  of  this  research.  The  calibration,  and  the  details  of  the 
research  program  and  data  for  this  period  of  research  are  given  next. 

2.2.1  CALIBRATION  OF  THE  WIEDEMANN  FATIGUE  RESEARCH  MACHINES 

Calibration  procedures  of  the  Wiedemann  fatigue  research 
machines  are  presented  next  in  two  parts;  Wiedemann  (Modified)  and 
Wiedemann  (Unmodified) . 

2. 2. 1.1  WIEDEMANN  (MODIFIED) 

An  analysis  was  performed  to  determine  the 
correct  loading  schedule  required  to  apply  a  desired  maximum, nominal 
bending  stress  to  specimens  of  this  research.  This  investigation 
required  that  freebody  diagrams  be  constructed  for  the  components 
of  the  system  using  actual  measured  weights  and  forces.  The 
components  of  the  system  consist  of  the  following: 

a)  The  auxilliary  upper-pan  counter-weight  system.  Fig.  2.2-1 

b)  The  bearing  and  housing  assembly.  Fig.  2.2-2 

c)  The  lower  pan  weight  system.  Fig.  2.2-3 

d)  Specimen  loading  freebody  diagram.  Fig.  2.2-4 

The  equations  of  force  were  written  for  each  of  the 
above  component  systems  treated  as  freebody  diagrams.  By  Figs.  2.2-1 
thru  2.2-4  ,  and  their  attendant  calculations  the  validity  of  the 

previous  calculations  provided  in  ONR  Report  2  [2,  pp.  67-83],  were 
ascertained . 


Fig. 2.2-1  Upper  Pan  Counter-Weight  System. 


Wp  =  pan  weight  =  0.067  pounds 

W  =  wire  weight:  W  =  0.0028  pounds,  W  =  0.0018  pounds 
c  Cj  c2 

W^  =  bracket  weight  =  0.068  pounds 

W*  =  pan  added  weight,  pounds 

When  unloaded  W'  =  0,  and  for  l  F  =  0 


W  +  If  )  -  (W,  •+  V  ) 


or 


0, 


(C.067  +  .0028)  -  (0.068  +  0.0018)  =  0, 
therefore 


(2.2-1) 


0.0698  -  0.0698  =  0. 

Thus  the  system  depicted  in  Fig.  2.2-1  is  in  balance  and  does  not  contribute 
to  the  specimen  loading  . 


I**-*  L . 

W/2 


W 


Fig.  2.2-2  Bearing  Housing  Assembly. 

length  of  moment  arm  for  lower  pan  load 

length  of  moment  arm  for  upper  pan  load 

diameter  of  specimen  (test  section) 

lower  pan  load 

upper  pan  load 

applied  stress  in  specimen 


V 


Fig.  2.2-3  Lower  Pan  Weight  Assembly. 


W„ 


weight  of  lower  pan  and  linkage,  5.960  pounds 


W, 


P  added  weight  to  lower  pan,  pounds 

r\ 


W/2 


load  applied  to  each  bearing  housi 


”g 


For  2  Fy  =  0,  (Fig. 2.2-3), 


2  (W/2)  -  (Wp  +  Wp  )  =  0 


or 


W  =  Wp  +  W_  . 

L  lA 


Since  Wp  was  determined  to  be  5.960  lb  , 
L 


w  =  5.960  +  W_ 


(2.2-2) 


(2.2-3) 


.1 


) 


T,  M 
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Li 


Fig.  2.2-4  Specimen  Loading  Free-body  Diagram. 

W’  =  upper  pan  added  weight,  pounds  - 

L’  =  upper  pan  moment  arm  =  4.65S  in. 

Wj,  =  weight  of  bearing  housing  at  C.  G.  =  4.219  lb.  | 

X  =  bearing  housing  moment  arm  =  2.074  in.  . 

Wp  =  weight  of  lower  pan  and  linkage  =  5.960  lb.  ! 

L 

L  =  lower  pan  and  linkage  moment  arm  =  4.000  in.  • 

Wp  =  lower  pan  added  weight,  pounds 

A  '  ! 

1  « 

W  =  weight  of  specimen  (1/2  specimen)  =  0.0358  lb. 

p  i 

Lq  =  specimen  moment  arm  =  4.967  in.  J 

bP 

; 

=  reaction  force  at  bearing,  pounds  j 

M  =  bending  moment  at  specimen  test  section,  in  -  lb. 

i 

I 


Determination  of  M,  the  Bending  Moment,  at 
the  Specimen  Groove 
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From  =0  (Fig.  2.2-4) 


CWp  +  Wp  ) 

rT  rA  ttft 

WHX  +  - - -  L  +  Ws  (Lg  )  -  L'  -  M  =  0, 


4.219(2.074)  + 


(5.960  +  W  ) 
*A 


4  +  0.0358(4.967)  -  —  (4.655)  =  M. 


M  =  8.750  +  11.920  +  2  W  +  0.178  -  2.328  W‘. 

A 


Therefore 


M  =  20.85  +  2  1  -  2.33  W’.*  * 

FA 


(2.2-4) 


Determination  of  the  Nominal  Stress  in  a  Specimen 


The  maximum  fiber  bending  stress  in  the  specimen  test  section  is  given 


(2.2-5) 


where 


I/C  = 


for  round  test  specimens  and  a  is  the  test  section  diameter.  Substitution 
for  I/C  in  Eq.  (2.2-5)  and  rearrangement  yields 


(2.2-6) 


Substitution  of  Eq. (2.2-4)  into  Eq.  (2.2-6)  yields 
20.85  +2  W  -  2.33  W 

s  = - li - , 

ud3/32 


*  All  figures  were  rounded  to  the  second  decimal  place. 
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1 

i 

i 


or 

212.45  +  20.38  Wp  -  23.74  W* 

s  =  - A -  .  (2.2-7) 

d-5 

Thus,  if  an  ungrooved  specimen  of  diameter  d  is  to  be  tested, 

pan  weights,  W  and  W,  can  be  chosen  to  give  the  desired  bending 
FA 

stress  using  Eg.  (2.2-7). 

To  calculate  the  upper  pan  weight,  W,  corresponding  to  various 
stress  levels  the  following  form  of  Eq.  (2.2-7)  is  useful: 

W  =  8.95  +  0.86  W  -  0.042  S  d^.  (2.2-8) 

FA 

In  order  to  reduce  the  number  of  variables  in  Eq.  (2.2-8)  and  to 
simplify  calculations,  the  lower  pan  added  weight  can  be  held  constant 
and  the  upper  pan  added  weight  varied  to  obtain  the  desired  bending 
stress.  The  lower  pan  added  weight  was  chosen  to  be  11.00  pounds; 
thus  Eq.  (2.2-8)  becomes 

W  =  8.95  +  0.86(11.00)  -  0.042  S  d3, 
or 

W  =  18.39  -  0.042  S  d3.  (2.2-9) 

Eq.  (2.2-9)  can  be  used  to  determine  the  upper  pan  weight  given  the 
bending  stress  and  the  specimen  test  diameter  for  ungrooved  specimens. 

Determination  of  the  Actual  Stress  in  a  Grooved  Specimen 

If  the  specimen  to  be  tested  is  grooved,  the  expression  which 
gives  the  pan  weight  is  needed.  The  stress  concentration  factor  must 
be  taken  into  account,  and  if  fatigue  is  involved  the  fatigue  notch 
factor  must  also  be  considered. 

The  procedure  is  as  follows: 

1.  Obtain  the  theoretical  stress  concentration  factor,  K^,  for 
the  notch  configuration  and  material  employed. 

2,  Obtain  the  fatigue  notch  factor  from 


-  —  *  „ - 


j 


(2.2-10) 


where  q,  the  notch  sensitivity,  and  Kt  are  obtained  from  published 
tables  and  figures  such  as  Peterson's  Stress  Concentration  Design 

fa££ors  [5;  p.  9,  Figs.  8  and  9;  pp.  48-49,  Figs.  39-40].  In  this 
report,  please  see  Figs.  2.2-5  thru  2.2-8. 


x  rjftj&ir-  I  ~  t~i-ri  \rr 


NOTCH  SENSITIVITY  FACTOR,  q 


FILLET  or*  HOLE  KO'JS.  f.  ACHES 


Fig.  2.2-5  Notch  sensitivity  of  normalized  steel  specimens 
[5,  p.  9] 


'NOTCH*  RADIUS.  •  "JCHES 

Fig.  2.2-6  Notch  sensitivity  of  quenched  and  tempered  alloy  steel 
specimens  (Based  on  Data  of  ^ough)  [5,  p.  9] 


ress  concentration  factor,  K^,  for  a  grooved  shaft  in  bending  (r/D  from  0.050  to  1.000) [5,  p.  49] 
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3.  Obtain  an  estimate  of  the  stress  in  the  groove  from 

S  =  K.  S  ,  (2.2-11) 

g  f  gn 

where  S  ,  the  nominal  grooved  stress,  is  calculated  by  using  Eq.  2.2-6. 

If  the  specimen  configuration  is  known  then  the  pan  loading  for 
a  desired  bending  stress  may  be  calculated.  Since  the  fatigue  notch  factor 
Kp,  is  greater  than  one,  the  actual  grooved  stress  will  be  greater  than 
the  nominal  grooved  stress.  It  should  be  clear  that  Eq. (2. 2- 10) gives 
only  an  estimate  of  the  fatigue  notch  factor  so  that  Eq.  (2.2-11)  will 
also  give  only  an  estimate  of  the  actual  grooved  stress.  To  find  the 
actual  K^,  two  Jots  of  specimens  are  tested:  One  lot  grooved  and  one 
lot  ungrooved;  however,  both  lots  have  the  same  groove  base  diameter. 


Then 


S 

nominal  for  ungrooved  specimens 

S 

nominal  for  grooved  specimens 


(2.2-12) 


In  Eq.  (2.2-12)  S  -  is  either  the  endurance  limit  or  endurance 

n  nominal 

strength  at  an  arbitrarily  large  number  of  cycles. 


Eq.  (2.2-12)  gives  the  actual  value  of  the  fatigue  notch  factor 
to  be  used  in  actual  stress  and  notch  sensitivity  calculations.  So 
when  testing  grooved  specimens  divide  the  nominal  ungrooved  stress  by 
the  estimated  fatigue  notch  factor  to  obtain  the  stress  to  be  used  in 
Eq.  (2.2-9). 

Thus 

W  =  18.39  -  0.042  S  d3.  (2.2-13) 

Kf 

Eq.  (2.2-13)  gives  the  desired  pan  weights  for  calculating  stresses 
for  grooved  specimens . 


M  L  j£.t. . 


r 

j 
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Thumb  Rules  for  Estimating  Starting  Stress  Levels 
for  Staircase  Testing 

The  method  of  staircase  type  testing  to  obtain  endurance 
limits  of  grooved  and  ungrooved  specimens  is  a  good  one.  However, 
the  time  involved  and  the  expense  of  manufacturing  test  specimens  are 
critical.  Therefore,  it  is  worthwhile  to  obtain  thumb  rules  for 
estimating  the  best  staircase  test  starting  point  so  that  realistic 
initial  stress  levels  are  employed  and  the  number  of  test  specimens 
used  is  minimized. 

After  investigating  staircase  test  data  for  ungrooved  specimens 
it  was  found  that  a  starting  stress  level  of  about  0.54  of  the  ultimate 
strength  is  the  optimum  starting  point.  Therefore,  when  staircase 
testing  ungrooved  specimens,  a  starting  stress  level  of 

S  .  =  0.54  SfM_  „  (2.2-14) 

st  Ultimate  v  3 

is  recommended  to  be  used  in  Eq.  (2.2-9)  in  place  of  S. 

For  staircase  testing  grooved  specimens  the  following  relation¬ 
ship  is  suggested  for  a  starting  stress: 


S  =  0.54  S„_  .  . 

st  Ultimate  . 


(2.2-15) 


This  stress  is  recommended  to  be  used  in  Eq.  (2.2-13)  in  place  of  S  . 

COMPUTER  PROGRAM  TO  DETERMINE  THE  WEIGHT  TO  BE  ADDED  TO  THE 
UPPER  PAN  TO  OBTAIN  THE  DESIRED  SPECIMEN  NOMINAL  STRESS 

A  computer  program  and  accompanying  output  is  given  in  Appendix 
A.*  This  output  gives  the  added  pan  weight  needed  to  achieve  a 
desired  stress  level  in  the  base  diameter  of  test  specimens.  Estimates 
of  the  groove  stress  and  pan  weights  to  be  added  are  also  given.  It 
is  recommended  that  this  program  be  consulted  before  any  testing  takes 
place. 


"Excerpts  from  Appendix  A  are  given  in  Table  2.2-1. 


,  j-  ,v--  ‘  J/ted^*^-^***'***' 


Table  2.2-1  Computer  Output  (Sample)  for  Pan  Weight  Calculations 
of  Wiedemann  Fatigue  Machine  (Modified)?' 


1.  Computer  program  can  be  found  in  Appendix  A. 

2.  The  groove  base  diameter,  d,  is  0.0937  in  for  all  groove  radii. 

3.  RAD  =  specimen  groove  radius. 

4.  CKF  =  dynamic  stress  concentration  factor. 

5.  SN  =  stress  in  base  of  groove. 

6.  SU  =  nominal  stress  in  base  of  groove. 

7.  WUPA  =  upper  pan  weight  to  be  added. 


Use  Appendix  A  as  follows: 
Ungrooved  Specimens 
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1.  For  ungrooved  specimens,  find  the  value  of  the 

material  and  calculate  S  using  Eq.  (2.2-14)  or  S  =  0.54 

^Ultimate' 

2.  Enter  the  Table  in  Appendix  A  whose  first  column  has  values 
RAD  =  1.870,  or  a  groove  radius  of  1.870  in.  which  results 
in  practically  no  stress  concentration,  hence  it  is  equi¬ 
valent  to  no  groove.  Go  down  the  SU  column  in  this  table 
and  stop  at  the  value  closest  to  that  found  in  Step  1,  then 
read  off  the  value  of  WUPA  in  the  adjacent  column  at  the 
right.  This  is  the  weight  that  has  to  be  added  to  the 
upper  pan,  with  the  lower  pan  weight  fixed  at  11.00  pounds, 
to  obtain  the  starting  stress  level  for  staircase  testing, 
and  is  the  value  given  by  Eq.  (2.2-9)  with  S  =  0.54 

S.,-..  *  substituted  in  it. 

Ultimate 

3.  The  SU  value  used  in  Step  2  to  arrive  at  the  NUPA  value  is 
the  stress  to  be  entered  in  the  staircase  plot. 

4.  If  the  specimen  does  not  fail  at  or  before  the  chosen  number 
of  cut-off  cycles  of  3  x  106  cycles,  run  a  fresh  specimen 
at  a  stress  level  2,000  psi  higher.  Do  this  by  getting  the 
1VUPA  v  lue  for  SU  2,000  psi  higher  than  that  used  in  the 
specimen  just  run.  If  the  specimen  fails  before  the  chosen 
number  of  cut-off  cycles,  run  a  fresh  specimen  at  a  stress 
level  2,000  psi  lower.  Do  this  by  getting  the  WUPA  value 
for  2,000  psi  lower  than  that  used  in  the  specimen  just  run. 
Continue  the  test  until  35  useful  specimens  are  run.  Have 
the  test  results  checked  by  your  supervisor  prior  to  term¬ 
inating  the  tests. 
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Grooved  Specimens 


1.  For  grooved  specimens  find  the  ultimate  value  of  the 
material  and  calculate  S  using  Eq.  (2.2-15). 

2.  Calculate  using  Eq.  (2.2-10)  after  getting  an  estimate 
for  q,  for  the  material  to  be  tested,  from  Figs.  2.2-5  and 
2.2-6  and  for  using  Figs.  2.2-7  and  2.2-8.  Such 
values  are  given  in  the  second  column  of  Table  2.2:1 
headed  CKF. 


3. 


4. 


Divide  the  S  value  found  in  Step  1  by  the  K_  value 

gn  y  1  f 

found  in  Step  2,  or  find  S  /K_. 

gn  f 

Enter  the  Table  in  Appendix  A  where  the  first  column  has 
values  RAD  equal  to  the  groove  radius  of  the  specimens  to 
be  tested.  Go  down  the  SU  column  and  stop  at  the  value 
closest  to  that  found  in  Step  3,  then  read  off  the  value 
WUPA  in  the  adjacent  column  at  the  right.  This  is  the 
weight  that  has  to  be  added  to  the  upper  pan  to  obtain  the 
starting  stress  for  staircase  testing.  This  value  is  the 
one  that  would  be  obtained  if  Eq.(2.2.-13)  were  used. 


5.  The  SU  value  used  in  Step  4  to  arrive  at  the  WUPA  value  is 
the  stress  to  be  entered  in  the  staircase  plot. 

6.  If  the  specimen  does  not  fail  at  or  before  the  chosen  number 
of  cut-off  cycles  of  3  x  10^  cycles,  run  a  fresh  specimen 

at  a  stress  level  2,000  psi  higher.  Do  this  by  getting  the 
NUPA  value  for  SU  2,000  psi  higher  than  that  used  in  the 
specimen  just  run.  If  the  specimen  fails  before  the  chosen 
number  of  cut-off  cycles,  run  a  fresh  specimen  at  a  stress 
level  2,000  psi  lower  than  that  used  in  the  specimen  just 
run.  Continue  the  test  until  35  useful  specimens  are  run. 
Have  the  test  results  checked  by  your  supervisor  prior  to 
terminating  the  tests. 
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CONCLUSIONS 

This  study  verifies  that  the  calibration  and  the  loading 
schedule  appearing  on  ONR  Technical  Report  4  [4,  pp.  212-220]  are 
correct . 
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2. 2. 1.2  WIEDEMANN  (UNMODIFIED) 

An  analysis  was  performed  to  determine  the  correct 
loading  schedule  required  to  apply  a  desired  maximum  bending  stress 
to  specimens  of  this  research  for  the  unmodified  machine. 

The  first  step  in  the  calibration  procedure  was 
to  weigh  the  bearing  housings,  loading  pan  and  pan  linkage  so  that 
the  contribution  of  each  component  to  the  bending  moment  in  the 
test  specimen  could  be  calculated. 

The  next  step  in  the  calibration  was  to  determine 
the  moments  due  to  the  bearing  housings  and  loading  linkage.  The 
procedure  was  to  support  the  bearing  at  its  normal  pivots  and  measure 
the  force,  F,  exerted  at  the  end  of  the  housing  shown  in  Fig.  2.2-9  • 
Knowing  the  distance  from  the  pivot  to  the  end  of  the  housing,  L, 
the  moment  in  the  specimen  was  calculated  for  the  housing  only. 

The  moment  due  to  the  linkage  and  load  pan  was  obtained  by  taking 
the  product  of  one-half  the  weight  of  the  linkage  and  the  distance 
from  the  knife  edges  to  the  pivot  point. 

Determination  of  Bending  Moment 
Due  to  the  Bearing  Housing 

The  bending  moment,  ,  at  the  specimen  groove  due  to  the  bearing 
housing  is 


(F  +  W  )  L, 
bP 


(2.2-16) 


1  Bq  =  bendin»  rrnmnt  due  to  the  hearing  housing  ami 
one-half  of  specimen  weight,  in.  -  lb 

F  -  force  component  due  to  bearing  housings  =  1.6  lb 
IVSp  =  WeiZht  of  cme-half  test  specimen  =  0.0358  lb 


L  -  moment  arm  for  force  =  4.655  in. 
Substitution  of  the  above  values  into  Eq.  (2.2-16)  gives 
»Bo  =  (1.6  +  0.0358)  (4.665), 

K  =  1.6358  X  4.655, 

0 


Determination  of  Sending  Moment 
Due  to  the  Pan  and  Linkage 

The  bending  moment  due  to  the  weight  of  the  pan  and 
accompanying  linkage  is  calculated  from 


(2.2-17 


where 

MPQ  =  bending  moment  at  the  specimen  due  to  the  pan  and 
linkage,  in.  -  lb 

P  =  weight  of  pan  +  linkage  =  1.575  +  4.375  =  5.950  lb 

Lp  “  length  of  moment  arm  for  pan  load  =  4.00  in. 


Substitution  of  these  val 


ues  into  Eq.  (2.2-17)  yields 
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Mp  =  11.90  in.  -  lb  . 

*0 

Thus  the  bending  moment  in  the  specimen  contributed  by  the 
pan  and  linkage  is  11.90  in.  -  lb. 

Determination  of  the  Total  Bending  Moment 

Due  to  the  Bearing  Housings,  Pan 

and  Linkage 

The  total  bending  moment  contributed  by  the  bearing  housings, 
the  pan,  and  the  linkage  is  calculated  as 

Mr  +  Mp  ,  (2.2-18) 

0  0 

the  total  bending  moment  due  to  the  bearing  housings 
and  the  pan-linkage  assembly,  in.  -  lb 

7.61  in.  -  lb 

11.90  in.  -  lb 

Substitution  of  these  values  into  Eq.  (2.2-18)  yields 

NL  =  7.61  +  11.90, 

0 

f-L,  =  19.51  in.  -  lb 

‘o 

Determination  of  the  Equivalent  Pan  Load 
The  equivalent  pan  load,  P  ,  may  be  defined  as  the  force  acting 

v 

at  the  knife  edge  supports  of  the  bearing  housings  due  to  the  weight 
of  the  bearing  housings  and  pan-linkage  assembly.  It  is  this  force 
that  produces  a  bending  moment,  M^,  in  the  test  specimen  prior  to 
loading  of  the  pan  and  testing.  It  is  necessary  to  determine  this 
load  so  that  the  actual  stress  produced  in  the  test  specimen  by  the 
addition  of  a  pan  load  may  be  determined.  The  equivalent  pan  load 
may  be  determined  as  follows: 


T 


0 


where 


o 


Mr 


0 


M 


e 


(2.2  -19) 


where 

i'L  =  19.51  in.  -  lb 
0 

Lp  =  4.00  in. 

Substitution  of  these  values  into  Eq.(2.2-19)  yields 

P  -  2C19.51) 

6  4 

Pe  =  9.76  lb. 

P.gg£rffination  of  the  Center  of  Gravity  of  the  Bearing  Housing 


The  center  of  gravity,  X,  may  be  determined  by  taking  the 
summation  of  moments  around  the  pivot  point  of  the  bearing  assembly 
as  follows: 


0, 


WD  CX) 
UB 


0, 


(2.2-20) 


where 

Wp  =  weight  of  the  bearing  housing  =  4.221  lb 
B 

X  =  center  of  gravity,  in. 

M-  =  7.61  in.  -  lb  . 

°0 


Substitution  of  these  values  into  Eq.  (2.2-20) 
7.61  -  4.221  (X)  =  0, 
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X  =  1.804  in. 

Determinntion  of  Stress  in  an  IJngrooved  Specimen 

Knowing  the  total  pan  weight  and  its  accompanying  moment  arm 
it  is  now  possible  to  calculate  the  stress  in  the  ungrooved  test 
specimen  from 


where 


I/C 

'  32  * 


for  round  test  specimens  and  d  is  the  specimen  test  section  diameter. 
Substitution  of  Eq.  (2.2-22)  into  Eq.  (2.2-21)  yields 


32  M 

3  * 

nd 


(2.2-25) 


Now  substitution  of  M  into  Eq.  (2.2-23)  gives 


S  - 


V  lp 

TTd3 


(2.2-24) 


where 

P^,  is  the  total  pan  weight  or  the  pan  added  weight  plus  the 
equivalent  pan  load.  Use  of  this  definition  in  Eq.  (2.2-24)  yields 

32 (-—)  (iv'p  +  9.76)  4.00 


Determination  of  Added  Pan  Weight  -  Stress  Relationship  for 
Un grooved  Specimen 

It  is  desired  to  establish  a  relationship  for  determining  the 
amount  of  weight  to  be  added  to  the  pan  to  achieve  the  desired 
maximum  bending  stress.  Rearrangement  of  Eq.  (2.2-25)  yields 


nd3S 

64 


9.76. 
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Thus  Eq.  (2.2-26)  provides  a  way  to  calculate  the  needed 
added  pan  weight  to  give  the  desired  maximum  bending  stress  in 
ungrooved  specimens. 

Determination  of  the  Maximum  Stress  in  a  Grooved  Specimen 

If  the  specimen  to  be  tested  is  grooved,  the  expression  which 
gives  the  pan  weight  for  a  desired  maximum  stress  in  the  specimen 
groove  is  needed.  The  stress  concentration  factor  must  be  taken 
into  account,  and  if  fatigue  is  involved  the  notch,  factor  must  also 
be  considered. 


The  procedure  is  as  follows: 

1.  Obtain  the  theoretical  stress  concentration  factor,  Kt> 
for  the  groove  configuration  and  type  of  loading  employed. 

2.  Obtain  the  fatigue  notch  factor  from 


Kf  =  q(Kt  -  1)  +  1,  (2.2-27) 

where  q,  the  notch  sensitivity  factor,  and  Kt>  are  obtained  from 
published  tables  and  figures  such  as  those  in  Peterson  (5  ;  p.9. 
Figs.  8  and  9  and  p.  10;  Fig.  10;  pp.  47-50,  Figs.  38  thru  41]. 
These  are  given  here  as  Figs.  2.2-5  thru  2.2-8  and  Figs.  2.2-10 
thru  2.2-12. 


3.  Obtain  an  estimate  of  the  stress  in  the  groove  from 


S  =  S  , 
8  f  *n 


(2.2-28) 


where  S  ,  the  nominal  groove  stress,  is  calculated  by  using 
Eq.  (2.2-23). 

If  the  specimen  configuration  is  known  then  the  pan  loading  for 
a  desired  bending  stressmaybe  calculated.  Since  the  fatigue  notch 
fa. 'tor,  K^,  is  greater  than  one,  the  actual  groove  stress  will  be 
grea^'r  that  the  nominal  groove  stress.  It  should  be  clear  that 


rta  a*  rnstammavfHUii 


' 

J 
1  1 

, 

1 

i 
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0  0.02  0.04  0.08  0.08  0.10  0.12  Oi4  0.16  ai8  020  0.22  024  026  028  030 

Pig.  2.2-11  Stress  concentration  factor,  K ,  for  a  grooved  shaft  in  bending 
(r/d  from  1.01  to  °°)[S,  p.  47]. 


concentration  factor,  K  ,  for  a  grooved  shaft  in  bending (r/d  from  1.005  to  «  and 
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Eq.  (2.2-27)  gives  only  an  estimate  of  the  fatigue  notch  factor  so  that 
Eq.  (2.2-28)  will  also  give  only  an  estimate  of  the  actual  groove  stress. 
To  find  the  actual  two  lots  of  specimens  are  tested:  one  lot 
grooved  and  one  lot  ungrooved;  however,  both  lots  have  the  same  groove 
diameter. 


Then 


K 


f 


nominal  for  ungrooved  specimen 

S 

nominal  for  grooved  specimen 


(2.2-29) 


In  Eq.  (2.2-29)  S  .  is  either  the  endurance  limit  or  the 
n  nominal 

endurance  strength  at  an  arbitrarily  large  number  of  cycles. 

Eq.  (2.2-29)  however,  is  restricted  to  the  case  of  completely  reversed 
bending.  Eq.  (2.2-29)  gives  the  actual  value  of  K^.  to  be  used  in 
stress  and  notch  sensitivity  calculations.  So  when  testing  grooved 
specimens  divide  the  nominal  ungrooved  stress  by  the  estimated  fatigue 
notch  factor  to  obtain  the  stress  to  be  used  in  Eq.  (2.2-26). 


Thus 


ird5 
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9.76  . 


(2.2-30) 


Eq.  (2.2-30)  gives  the  desired  added  pan  weight  for  calculating  the 
stress  for  grooved  specimens. 


THUMB  RULES  FOR  ESTIMATING  STARTING  STRESS  LEVEL  FOR 
STAIRCASE-TESTING 


The  method  of  staircase  type  testing  to  obtain  endurance  limits 
of  grooved  and  ungrooved  specimens  is  a  good  one.  However,  the  time 
involved  and  the  expense  of  manufacturing  test  specimens  are  critical. 
Therefore,  it  is  worthwhile  to  obtain  thumb  rules  for  estimating  the 
best  staircase  test  starting  point  so  that  realistic  initial  stress 
levels  are  employed  and  the  number  of  test  specimens  used  is  minimized. 

After  investigating  staircase  test  data  for  ungrooved  specimens 
it  was  found  that  a  starting  stress  level  of  about  0.54  of  the  ultimate 


strength  is  the  optimum  starting  point.  Therefore,  when  staircase¬ 
testing  ungrooved  specimens,  a  starting  stress  level  of 
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S  =  0.54  SM1„.  +  , 

Ultimate 

is  recommended  for  use  in  Eq.  (2.2-26). 


(2.2-31) 


For  staircase-testing  grooved  specimens,  a  starting  stress 


level  of 


S  =  0.54  SM1  .  , 

g  Ultimate 


(2.2-32) 


is  recommended  for  use  in  Eq.  (2.2-30). 


Computer  Program  to  Determine  Pan  Added  Weights 


In  Appendix  B,  Table  2.2-2  will  be  found  a  computer  program 
and  accompanying  output.  This  output  gives  the  pan  weight  needed  to 
achieve  a  desired  stress  level .  Estimates  of  the  grooved  stress  and 
pan  weights  are  also  given.  It  is  recommended  that  this  program  be 
consulted  before  any  testing  takes  place.  Table  2.2-2  provides 
excerpts  from  Appendix  B. 


Use  Appendix  B  as  follows: 
Ungrooved  Specimens 


1. 


2. 


For  ungrooved  specimens,  find  the  Syjt^mate 
material  and  calculate  S  using  Eq.  (2.2-31), 


value  of  the 
or  S  =  0.54 


Q 

Ultimate' 

Enter  Appendix  B  table  whose  first  column  has  values 
RAD  =  1.870.  or  a  groove  radius  of  1.870  in.  which  results 
in  practically  no  stress  concentration,  hence  it  is  equi¬ 
valent  to  no  groove.  Go  down  the  SU  column  in  this  table 
and  stop  at  the  value  closest  to  that  found  in  Step  1,  then 
read  off  the  value  of  WPA  in  the  adjacent  column  at  the 
right.  This  is  the  weight  that  has  to  be  added  to  the  pan 


Table  2.2-2  Computer  Output  (Sample)  for  Pan  Weight  Calculations  of 
Wiedemann  Fatigue  Machine  (Unmodified).* 


RAD* ** 

CKF 

SN 

su 

WPA 

inches 

pounds 

.031 

1.510 

0 

0 

-9.76 

.031 

1.510 

1510 

-8.79 

.031 

1.510 

3020 

-7.83 

.031 

1.510 

4530 

■1 

-6.86 

.062 

0 

0 

-9.76 

.062 

1.390 

1390 

-8.79 

.062 

1.390 

2780 

-7.85 

.062 

• 

• 

1.390 

4170 

-6.86 

« 

• 

w&wm 

1.243 

0 

0 

-9.76 

Hill 

1.243 

1243 

-8.79 

.125 

1.243 

2486 

-7.83 

.125 

1.243 

3729 

• 

• 

H 

-6.86 

.250 

1.142 

0 

0 

-9  76 

.250 

1.142 

1142 

1000 

-8.79 

.250 

1.142 

2284 

2000 

-7.83 

.250 

1.142 

3426 

3000 

-6.86 

1  870 

1.000 

0 

0 

-9.76 

1.870 

1.000 

1000 

1000 

-8.79 

1.870 

1.000 

2000 

2000 

-7.83 

1.870 

1.000 

3000 

3000 

-6.86 

• 

• 

• 

• 

• 

' 

' 

*Complete  tables  are  found  in  Tables  10  .2-6  thru  10  .2-10  computer 
program  can  bo  found  in  Appendix  B. 

**The  groove  base  diameter,  d,  is  0.2700  in,  for  all  groove  radii. 
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to  obtain  the  starting  stress  level  for  staircase  testing, 
and  is  the  value  given  by  Eq.  (2.2-26)  with  S  =  0.54 

S  .  ^  substituted  in  it. 

Ulitraate 

3.  The  SU  value  used  in  Step  2  to  arrive  at  the  WPA  value  is 
the  stress  to  be  entered  in  the  staircase  plot. 

4.  If  the  specimen  does  not  fail  at  or  before  the  chosen 
number  of  cutoff  cycles  of  3  x  10^  cycles,  run  a  fresh 
specimen  at  a  stress  level  2,000  psi  higher.  Do  this  by 
getting  the  WPA  value  for  SU  2,000  psi  higher  than  that 
used  in  the  specimen  just  run.  If  the  specimen  fails  before 
the  chosen  number  of  cutoff  cycles,  run  a  fresh  specimen  at 
a  stress  level  2,000  psi  lower.  Do  this  by  getting  the  WPA 
value  for  2,000  psi  lower  than  that  used  in  the  specimen  Just 
run.  Continue  the  test  until  35  useful  specimens  are  run. 


Grooved  Specimens 


1.  For  grooved  specimens  find  that  ultimate  value  of  the  material 

and  calculate  S  using  Eq.  (2.2-32). 

^n 

2.  Calculate  K^.  using  Eq.  (2.2-27)  after  getting  an  estimate 
for  q,  for  the  material  to  be  tested,  from  Figs.  2.2-5, 

2.2- 6,  or  2.2-10,  and  for  using  Figs.  2.2-7,  2.2-8, 

2.2- 11,  or  2.2-12.  Such  values  are  given  in  the  second 
column  of  Table  2.2-2  headed  CKF. 


3. 


4. 


Divide  the  value  S  found  in  Step  1  by  the  K,.  value  found 
g„  f 

in  Step  2,  or  find  S  /K^. 

Enter  Appendix  B  table  where  the  first  column  has  values 
RAD  equal  to  the  groove  radius  of  the  specimens  to  be  tested. 
Go  down  the  SU  column  and  stop  at  the  value  closest  to  that 
found  in  Step  3,  then  read  off  the  value  WPA  in  the  adjacent 
column  at  the  right.  This  is  the  weight  that  has  to  bo 
added  to  the  pan  to  obtain  the  starting  stress  for  staircase 
testing.  This  value  is  the  one  that  would  be  obtained  if 
Eq.  (2.2-30)  were  used. 
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5.  The  SU  value  used  in  Step  3  to  arrive  at  the  NPA  value  is 
the  stress  to  be  entered  in  the  staircase  plot. 

6.  If  the  specimen  does  not  fail  at  or  before  the  chosen 
number  of  cutoff  cycles  of  3  x  10  cycles,  run  a  fresh 
specimen  at  a  stress  level  2,000  psi  higher.  Do  this  by 
getting  the  WPA  value  for  SU  2,000  psi  higher  than  that  used 
in  the  specimen  just  run.  If  the  specimen  fails  before  the 
chosen  number  of  cutoff  cycles,  run  a  fresh  specimen  at  a 
stress  level  2,000  psi  lower.  Do  this  by  getting  the  WPA 
value  for  2,000  lower  than  that  used  in  the  specimen  just 
run.  Continue  th<»  test  until  35  useful  specimens  are  run. 


CONCLUSIONS 

This  study  indicates  that  this  calibration  report  and  the 
associated  loading  schedule  should  replace  those  appearing  in  ONR 
Technical  Report  No.  2.  [2,  pp.  67-73]. 

RECOMMENDATIONS 


For  any  research  using  the  unmodified  R.  R.  Moore  Reliability 
Research  Machine  the  following  equations  are  recommended: 

1.  To  obtain  the  ungrooved  specimen  maximum  bending  stress  use 

32(~)0Vp  +9.76)  4.00 

S  =  - r -  (2.2-25) 


2.  To  obtain  the  added  pan  weight  for  a  desired  maximum  bending 
for  ungrooved  specimens  use 


7Td3S 


64 


9.76  . 


(2.2-26) 


I 


6S 


3.  To  obtain  the  added  pan  weight  for  a  desired  maximum 
bending  stress  for  grooved  specimens  use 


md^  S 

*» 

64  K 


9.76 


(2.2-30) 


4.  To  obtain  a  starting  stress  level  for  staircase  testing 
of  ungrooved  specimens  use  the  instructions  on  pages  71-74. 


5.  To  obtain  a  starting  stress  level  for  staircase  testing  of 
grooved  specimens  use  the  instructions  on  pages  74-75. 
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2.2.2  RESEARCH  PROGRAM  AND  DATA 

Tables  2. 5. -3  and  2.2r4  provide  a  cross-reference  for 
the  Wiedemann  machines  test  program,  and  the  lata  generated  and 
reduced.  The  test  program  in  these  tables  are  broken  down  by  data 
type,  i.e.,  Cycles-to-failure  and  Staircase  testing  data. 
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Table  2.2-3  Cyclcs-To-Fai lure  on  Wiedemann  Machines. 


Group 

No. 

Notch  Radius 

Inches  +  .001 
-  .002 

Sample 

Size 

Volume  II 

Section  Cross  Reference 

Raw  Data 
Table  No. 

Reduced  Data 

Table  No.  Figure  No. 

AISI  4130  Steel  Rod  D 

=  0.375 

in.  diameter 

d 

=  0. 

0937  in.  dia.  of  Notch 

82 

0.031 

35 

Test 

to 

be 

re-run 

83 

0.031 

35 

Test 

to 

be 

re-run 

84 

0.031 

35 

Test 

to 

be 

re-run 

85 

0.031 

35 

Test 

to 

be 

re-run 

104 

0.062 

35 

Test 

to 

be 

re-run 

105 

0.062 

35 

Test 

to 

be 

re-run 

106 

0.062 

35 

Test 

to 

be 

l'e-run 

107 

0.062 

35 

Test 

to 

be 

re-run 

88 

0.125 

55 

Test 

to 

be 

re-run 

89 

0.125 

35 

8.2-11 

9.2-1  to  9.2-3 

yo 

0.125 

35 

8.2-14 

9.2-4  to  9.2-6 

91 

0.125 

35 

8.2-16 

9.2-7  to  9.2-9 

92 

0.250 

35 

Test 

to 

be 

re-run 

93 

0.250 

35 

8.2-18 

9.2-10  to  9.2-12 

94 

0.250 

35 

8.2-20 

9.2-13  to  9.2-1S 
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Table  2.2-4  Notch  Sensitivity  Fatigue  Test  Data  -  Staircase  Testing 
on  Wiedemann  Machines 


Group 

No. 

Notch  Radius 
Inches  +  „00| 

Sample 

Size 

Volume  II 

Section  Cross 

Reference 

Raw  Data 
Table  No. 

Reduced  Data 

Table  No.  Figure  No. 

AISI  1038  Steel  Rod  D 

=  0.375  in.  diameter. 

d  =  0.2700  in.  dia.  of  Notch 

162 

0.031 

35 

8.2-23 

9.2-1 

9.2-19 

163 

0.062 

35 

8.2-25  ; 

9.2-2 

9.2-20 

164 

0.125 

35 

8.2-27 

9.2-3 

9.2-21 

165 

0.250 

35 

8.2-29 

9.2-4 

9.2-22 

166 

1.870 

35 

8.2-31 

9.2-5 

9.2-23 

2.3 


AXIAL  FATIGUE  RELIABILITY  RESEARCH  MACHINE 


2,3.1  AXIAL  FATIGUE  TESTING  TO  GENERATE  GOODMAN  DIAGRAMS 

The  axial  fatigue  machine  can  apply  a  constant  mean 
tensile  load  onto  which  a  constant  amplitude  alternating  tensile- 
compressive  load  can  be  superimposed  along  the  longitudinal  axis  of  a 
test  specimen.  The  application  of  the  combined  loads 
provides  the  capability  to  determine  the  endurance  strength  of  the 
test  specimens  an!  to  investigate  the  effects  of  interaction  between 
alternating  and  mean  loads.  When  the  endurance  strength  calculated  at 
specified  ratios  of  alternating  to  mean  stresses  are  plotted  on  a 
biaxial  graph,  a  strength  surface  known  as  the  distributional  Goodman 
diagram  results,  as  shown  in  Fig.  2.3-1. 

Four  different  theories  are  presented  in  Fig.  2. 3.-1. 

The  first  one  is  the  modified  Goodman  line  joining  the  mean  of  the 
endurance  strength,  at  a  specific  cycles  of  life,  to  the  mean  of  the 
static  ultimate  strength  of  the  material.  This  a  conservative  design 
criterion  for  combined-stress  fatigue  and  may  be  used  when  fracture  is 
the  failure  mode.  The  second  is  a  modification  of  the  previous  one, 
in  that  the  line  joins  the  means  of  the  endurance  strength  to  the  mean 
of  the  static  yield  strength  instead  of  the  ultimate.  This  is  also  a 
conservative  combined-stress  design  crite’ aon  and  may  be  used  i-hen 
yield  is  the  failure  mode.  The  third  criterion  is  the  Gerbe*.  parabola 
passing  through  the  mean  of  the  endurance  strength  and  the  mean  of  the 
static  ultimate  strength.  This  parabola  results  from  the  maximum 
principal  shear  stress  theory  of  failure  [6,  p.  183],  This  theory  is 
optimistic  relative  to  the  first  two  criteria  and  yet  is  closer  to  the 
actual  fatigue  strength.  This  criterion  is  not  recommended  for  stresses 
beyond  the  yield  strength.  The  fourth  criterion  is  an  ellipse, 
resulting  from  the  von  Mises-Hencky  failure  criterion  [7,  p.  85],  which 
passes  through  the  mean  of  the  endurance  strength  and  the  mean  of  the 
static  ultimate  strength.  This  criterion  is  the  closest  to  the  actual 
combined-stress  fatigue  strength  of  most  steels.  Again  this  criterion 
is  not  recommended  for  stresses  beyond  the  yield  strength,  however,  for 
lack  of  better  criteria,  it  is  being  used  to  represent  fracture  as  well. 


Alternating  stress 
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Fig.  2.3-1  Goodman  diagrams  showing  various  theories  defining  the 
safe  design  region  for  combined-stress  fatigue. 
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The  axial  fatigue  machine  was  used  to  obtain  experi¬ 
mental  data  to  determine  the  strength  distributions  for  a  specific 
cycles  of  life  for  specimens  of  a  specified  steel  and  geometry  for 
construction  of  distributional  Goodman  diagrams.  In  order  to  obtain 
sufficient  data  points  to  adequately  define  the  strength  surface,  axial 
fatigue  tests  are  conducted  at  stress  ratios  of  infinity,  2.0,  1.0,  0.4 
and  0.2.  The  static  ultimate  strength  distribution  is  used  for  the 
strength  distribution  at  a  stress  ratio  of  zero. 


2. 3. 1.1  DETERMINATION  OF  ENDURANCE  STRENGTH 

DISTRIBUTION  AT  THE  STRESS  RATIO  r,  =  ® 

_  ;  _  3 

First  an  estimate  of  the  endurance  strength 

mean  is  required.  Table  2.3.1  was  developed  to  obtain  the 

ratio  between  the  endurance  strength  mean  and  the  ultimate 

strength  mean,  based  on  experimental  data  for  the  four  materials 

researched.  The  results  indicate  that  the  value  S  /S  =0.45 

e  u 

is  a  representative  one.  Consequently  the  following  equation 
was  developed  for  the  estimate  of  the  mean  of  the  endurance 
strength  from  a  knowledge  of  the  mean  of  the  static  ultimate 
strength: 

S  =  0.45  S  .  (2.3-1) 

e  u  v 


Knowing  Sufrom  the  experimental  data  for  static  ultimate 
strength,  Sq  can  be  calculated  from  Eq.  (2.3-1).  The  starting 
load  for  the  staircase  test  may  then  be  obtained  from 


where 


(2.3-2) 


=  mean  alternating  load  in  lb. 

A  =  mean  cross  sectional  area  of  the  test  specimen 

.  .  2 
m  m  . 


Table  2.3-1  Ratio  of  experimental  endurance  strength  to  static 
ultimate  strength  for  staircase-test  starting-stress 
determination. 


Material 

.  Mean  of 
Static  Ultimate 
Strength,  S 

PSi  U 

Mean  of 
Endurance 
Strength,  S 
psi  c 

S  /S 
e  u 

AISI  1018 

60,000 

26,000 

0.434 

AISI  1038 

69,000 

31,000 

0.449 

AISI  4130 

105,000 

40,000 

0.381 

AISI  4340 

116,000 

49,000 

0.423 
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The  staircase  method  of  testing  requires  that  if  the  first 
specimens  fails  before  2  x  20°  cycles,  the  i.ext  specimen  is 
loaded  to  a  stress  level  one  increment  higher.  This  increment 
has  been  chosen  to  be  approximately  1,000  psi. 


2.3. 1.2  DETERMINATION  OF  STRENGTH  DISTRIBUTION  AT 
STRESS  RATIOS  OTHER  TRAN  « 


For  stress  ratios  other  than  »  it  is  necessary 
to  apply  a  mean  load  as  well  as  an  alternating  load.  This  is 
done  by  determining  the  mean  as  well  as  the  alternating  loads 
as  follows:, 

Previous  experimental  research  at  The  University 
of  Arizona  provides  substantial  evidence  that  the  actual 
strength  distribution  for  axial  fatigue  lie  in  between  the 
Gerber  parabola  and  the  von  Mises-Hencky  ellipse.  The 
equation  of  the  Gerber  parabola  is 


(2.3-3) 


and  the  equation  of  the  von  Mises-Hencky  ellipse  is 

(-TT-)  *  (r-)  -  »•  P-3-'') 

v  e  '  u 

The  data  appears  to  fa /or  the  ellipse,  furthermore  a  slightly 

higher  stress  vector,  Svm  value  will  provide  a  good  starting 

stress  with  the  first  few  specimens  being  failures.  Thus  as 

an  estimate  of  the  mean  of  for  starting  staircase  testing 

a  value  favoring  the  von  Mises-Hencky  ellipse  is  used  for  a 

2 

specific  stress  ratio.  Dividing  Eq.  (2.3-4)  by  (S  )  gives 


9 


+ 


(2.3-5) 


parabola 

=  1.0 


v°n  Mises-Hencky  ellips 


stress  Sector lrqC3^  ^  sh°win^  the  e  f 

'  V  at  diffe^t  sSsTS,-lf  the 
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S 

and  solving  for  S  and  substituting  r  =  ■= —  gives 
m  So 

m 


s 


1/2  ' 


(2.3-6) 


Eq.  (2.3-6)  yields  an  estimate  of  in  terms  of  known 
quantities. 


Now  stress  vector  S  may  be  calculated  from 


S  =  (S  2  +  S  2)i/2 
v  v  a  m  ' 


(2.3-7) 


where 


Alternating  load _ 

Cross  sectional  area  of  specimen 


> 


and 


S  * 
m 


Mean  load 


Cross  sectional  area  of  specimen 
S 


Substitution  of  r 


in  Eq.  (2.  /-7)  j  lelds 


m 


Sv  *  Sm  Crs2  *  1>1/2‘  P-3-*) 

Thus  the  value  of  S  found  from  Eq.  t.2.3-6)  and  applied  to  Eq.  (2,3-8) 

provides  an  estimate  of  the  mean  of  S  ,  at  a  specified  stress  ratio. 

Points  on  the  von  Mises-Hencky  ellipse  calculated  by  this  method1  and 

used  as  initial  values  of  S  for  staircase  test.'  at  various  stress 

v 

ratios  are  shown  in  Fig.  2.3. -2. 

The  estimates  of  S  and  S  must  be  translated  into  axial  mean 
m  v 

and  alternating  loads  for  testing  on  the  Axial  Fatigue  Research  Machine. 
After  Eq.  (2.3-6)  is  used  to  obtain  an  estimate  of  S  ,  the  required 
mean  load  is  determined  from 


P  =  .s 
m  m 


•  A  . 


(2.3-9) 


Footnote: 

1.  The  PDP-8  program  to  calculate  points  on  the  Gerber  parabola  and 
the  von  Mises-Hencky  ellipse  is  listed  in  Appendix  F. 
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i 

I 

i  - 


The  corresponding  alternating  load  is  then  calculated  from 


P 

a 


(2.3-10) 


The  mean  load  indicator  of  the  Axial  Fatigue  Machine  is  scaled 

in  one  pound  increments;  therefore,  the  estimated  values  of  S  ,  and 

m 

consequently  of  Sv,  are  adjusted  to  the  nearest  pound  of  force.  The 
revised  stress  values,  and  S^,  are  calculated  by  reverse  application 
of  Eqs .  (2.3-9)  and  (2.3-8). 


Staircase  testing  requires  steps  of  equal  increments  of  stress . 
For  this  program  an  acceptable  range  of  increments  of  Sv  is  1,000  to 
1,S00  psi.  Similar  to  revising  the  initial  stress  level,  a  method  is 
required  to  adjust  the  initially  selected  stress  increment  to  one 
which  corresponds  to  an  increment  in  multiples  of  pounds  of  mean  load. 
The  methodology  used  for  calculating  the  starting  stress  level,  S^, 
and  the  stress  increment,  is  given  below: 

(1)  Determine  an  estimate  of  S^  from  the  static  ultimate 
strength  of  the  specimen,  Eqs.  (2.3-1)  and  (2.3-6)  and  the 
desired  stress  ratio. 

(2)  Calculate  an  initial  value  of  Sv  from  Eq.  (2.3-7). 

(3)  Find  the  required  mean  load,  P^,  from  Eq. (2.3-9). 

(4)  Round  to  the  nearest  pound  and  revise  the  estimated 

value  of  S  . 

m 

(5)  Calculate  the  corresponding  values  of  S  and  S  for  the 

V  3. 

desired  stress  ratio. 


(6)  Select  a  tentative  value  for  the  stress  increment,  Is> 

(7)  Calculate  the  comparable  increment  in  mean  stress  from 
Eq.  (2.3-8)  rearranged  in  the  form 


I  =  I  /(r  2  +  1)1/2. 
m  s  ^  s  ' 


(2.3-11) 


(8)  Find  the  required  increment  in  mean  load,  I  ,  from  an 

*  m 

adaptation  of  Eq.  (2.3-9)  as 
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1 


(2.3-12) 


(9)  Round  the  calculated  value  of  I  to  the  nearest  pound, 

and  apply  Eqs.  (2.3-11)  and  (2.3-12)  in  reverse  order  to 
determine  the  reversed  stress  increment. 

(10)  Calculate  the  required  meair  and  alternating  loads  for  each 
level  of  anticipated  for  the  test  from  the  above 
equations.  Each  stress  level  of  the  test  is  determined 
from  the  preceding  stress  level  (depending  upon  whether 

g 

the  specimen  successfully  completed  2  x  10  cycles  or 
failed  prior  to  that  time)  by  the  equation 


=  S  ±  I 
v  s 


(2.3-13) 


A  PDP-8  computer  program  has  been  deveToped  to  calculate 
the  initial  stress  level  and  subsequent  stress  levels, 
and  the  required  mean  and  alternating  loads.  A  copy  of 
the  program  for  the  initial  test  point  is  listed  in 
Appendix  G,  and  for  subsequent  points  in  Appendices  H  and  I. 

2.3.!.5  PLOTTING  AND  REDUCTION  OF  STAIRCASE  DATA 

The  data  obtained  from  staircase  testing  are 
recorded  and  reduced  using  the  methods  discussed  for  the 
Wiedemann  machines  in  Section  3.2.3. 

2.3.2  RESEARCH  PROGRAM  AND  DATA 

During  this  reporting  period  AISI  1018  steel  specimens 
were  tested  to  get  the  strength  distributions  in  terms  of  the  mean 
and  the  standard  deviation  of  S^  at  different  stress  ratios.  Table 
2.3-2  summarizes  the  test  program  and  points  out  the  sections  of  this 
report  where  the  raw  data  and  the  reduced  data  can  be  obtained. 


Table  2.3-2  Axial  Fatigue  Machine  test  program  for  AISI  1018  steel 
specimens  and  references  for  the  raw  and  the  reduced 
data  for  groups  of  35  specimens. 


Group 

No. 

Stress 

Ratio 

Diameter  of 
test  section 
inches 

Raw 

Data 

Reduced  Data 
Table  No .  Figure 

161 

CO 

0.075 

8.3-5 

9.3-3 

9.3-3 

160 

2.0 

0.075 

8.3-4 

P.3-2 

9.5-2 

159 

1.0 

0.075 

8.3-3 

9.3-1 

9.3-1 

158 

0.4 

0.070 

Testing 

incomplete 

157 

0.2 

0.070 

Testing 

incomplete 

156 

0.0 

0.250 

8.3-1 

Not  applicable 
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3.0  DATA  REDUCTION  AND  RESULTS 

Data  generated  by  the  experimental  research  discussed  in  Chapter  2 
were  reduced  and  analyzed  to  determine  the  dynamic  strength  characteristics 
of  specimens  of  different  AISI  steels  and  geometry.  Cycles-to-failure  tests 
conducted  at  specified  stress  levels  were  analyzed  to  determine  which 
distribution  model  provided  the  best  fit  to  the  data  and  to  estimate  the 
parameters  of  the  distributions  being  evaluated.  Stress-to-failure  test 
data  were  reduced  to  determine  the  distribution  parameters  of  endurance 
strength  for  the  normal  distribution. 

3.1  CYCLES-TO-FAILURE  DATA 

The  cycles-to-failure  data  presented  in  Section  II  were  reduced  with 
the  aid  of  two  CDC-6400  computer  programs,  program  CYTOFR  and  program 
WEIBULL.  The  primary  functions  of  theprograms  are  to  calculate  distribution 
parameters  from  the  data  and  to  perform  goodness -of-fit  tests  to  provide 
indicators  of  the  best  distribution  model. 

3.1.1  PROGRAM  CYTOFR 

Program  CYTOFR  was  discussed  and  the  FORTRAN  program  statements 
were  presented  in  the  Fourth  ONR  Report  (  4  ,  pp.  163-174).  During  the 
past  year  a  number  of  changes  have  been  made  to  improve  its  effectiveness. 
The  current  program  is  listed  in  Appendix  C.  Program  CYTOFR  determines 
parameter  estimates  of  the  following  distributions: 


Normal : 

fNm 


a  /2tt 


having  the  parameters 

T  =  mean  of  cycles  to  failure 
o  =  standard  deviation  of  cycles  to  failure 
T  =  cycles  to  failure 


(3.1-1) 


Lognormal: 
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(T) 


1 

TV  /I7 

o 


T  -  T* 


(3.1-2) 


T'  =  mean  of  logg  T 

o'  =  standard  deviation  of  loge  T 

T  =  cycles-to-failure. 

Program  CYTOFR  performs  additional  calculations  and  makes  goodness- 
of-fit  tests  as  follows: 


1.  Coefficient  of  skewness,  a_. 

2.  Coefficient  of  Kurtosis,  • 

3.  Chi-Squared  Goodness-of-Fit  test. 

4.  Kolmogorov-Smirnov  Goodness-of-Fit  test. 

The  calculation  of  the  distribution  parameters  are  accomplished 
within  the  computer  program  by  incorporation  of  the  following  equations: 

The  data  mean,  T,  is  calculated  from 


n 


n 

where  T.,  i  =  1,  2,  — ,n  are  the  values  for  n  data  points. 

2 

The  unbiased  estimate  of  the  distribution  variance,  a  ,  is 
calculated  from 


(3.1-3) 


2 

a 


n  -  2 
l!lCTi  -  T) 

n-1 


(3.1-4) 


The  coefficients  of  skewness  and  kurtosis,  as  discussed  by  Hahn  and 
Shapiro,  are  summarized  by  the  moments  of  the  distribution  about  the  mean 
connected  to  dimensionless  coefficients  for  simplicity  of  use  [8  ,  pp. 
45-49].  The  coefficient  of  skewness,  abuses  the  relationship  between 
the  second  moment  which  measures  dispersion  and  the  third  moment  which 
measures  skewness  to  obtain  a  standardized  measure  of  skewness  of  the 
distribution  relative  to  its  degree  of  spread.  If  the  distribution  is 
normal  =0. 
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The  quantitive  measure  is  found  from  the  equation 


where 


and 


(5.1-5) 

(3.1-6) 

(3.1-7) 


Similarly  the  coefficient  of  kurtosis,  ci^,  uses  the  relationship 
between  the  second  moment  and  the  third  moment  of  the  distribution  about 
the  mean  to  obtain  a  standardized  measure  of  peakedness  of  the  distri¬ 
bution  relative  to  its  degree  of  spread.  The  quantitative  value  of  a^, 
which  is  equal  to  3.0  for  an  exact  normal  distribution,  is  obtained  from 
the  equation 


where 


(3.1-8) 


(3.1-9) 


and 

M2  is  as  given  in  (3.1-7  ). 

The  Chi-Squared  goodness-of-fit  test  is  the  comparison  of  actual 
observations  of  grouped, ordered  data  with  the  expected  number  of  observations 
calculated  using  the  parameter  estimates  and  the  appropriate  distribution 
equation.  The  methodology  for  application  of  this  test  is  contained  in  the 
Third  Technical  Report  [  3  ,  pp.  256-267]. 
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The  Kolmogorov-Srairnov  goodness-of-fit  test  is  the  comparison 
of  the  cumulative  probability  density  graph  of  the  data  distribution 
calculated  from  the  estimated  parameters  with  the  relative  cumulative 
graph  of  each  ordered  data  point.  Discussions  of  the  application  of 
the  K-S  test  to  normal  and  lognormal  distributions  are  presented  in 
the  Fourth  Technical  Report  [4,  pp.  51-58]. 

The  normal  and  lognormal  distributions  and  their  applications 
to  design-by-rel 5 ability  methodology  for  dynamic  mechanical  components 
are  discussea  in  the  First  ONR  Report  [1,  pp.  44-5?]. 


3.1.2  PROGRAM  WE1BULL 


Program  WEIBULL  has  extended  the  data  reduction 
techniques  by  the  computerized  calculation  of  Weibull  distribution 
parameters  from  the  data,  and  the  application  of  the  X^'  and  K-S  goodness- 
of-fit  tests  indicated  above.  The  Weibull  distribution  is  described 
by  the  function 


fNm 


6-1 

e 


(3.1-10) 


where 

8  =  shape  parameter, 
n  =  scale  parameter, 

Y  =  location  parameter. 


Consideration  of  the  Weibull  distribution  as  an  appropriate 
model  for  the  distribution  of  fatigue  data  was  introduced  and  discussed 
in  the  Fourth  Technical  Report  [4,  pp.  63-84].  Development  of  the 
analysis  capability  was  limited  to  manual  techniques  until  program 
WEIBULL  was  completed. 

The  basic  FORTRAN  computer  program  to  solve  for  estimates  of 
Weibull  parameters  and  cycle  life  for  Specified  levels  of  reliability 
was  provided  by  Mr.  Thomas  C.  Stansberry,  Delco  Radio  Division,  General 
Motors  Corporation.  His  program  was  adapted  to  The  University  of 
Arizona  CDC-6400  computer  and  updated  to  include  subroutines  for  the 
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Chi -Squared  and  Kolraogorov-Smimov  goodness-of-fit  tests.  A  copy  of 
the  program  FORTRAN  statements  is  contained  in  Appendix  D. 

The  first  data  card  contains  the  sample  size  and  the  minimum 
life  increment  for  use  in  linearizing  the  x-y  relationship.  Subsequent 
data  cards  (one  for  each  specimen)  contain  the  failure  information. 

The  first  operation  performed  by  the  computer  is  to  establish  an 
ordered  array  of  the  failures  and  median  ranks  similar  to  standard 
tables  [4,  p.  67].  However,  the  computer  calculates  the  median  ranks 
so  that  manual  use  of  tables  for  each  data  sample  is  unnecessary.  _ 
The  next  operation  is  the  calculation  of  y.  =  In  In  - — ?- 

1  L1  “  F(ti3_ 

and  x^  =  ln(t^  -  y^) .  Where  i  =  1,  2,  — ,  n  and  y^  =  minimum  life 
increment  (1,  2,  — ,  k)  such  that  y.  <  t. .  As  the  array  of  y.  and 

K  X  X 

x,  is  computed  the  method  of  least  squares  is  used  to  determine  the 
degree  of  linearity.  This  operation  is  iterated  with  y^  being 
increased  in  increments  until  the  best  fit  straight  line  is  obtained. 

At  that  time  the  computer  records  the  estimates  of  8,  y,  and  n-  It 
then  calculates  the  one  percent  life,  ten  percent  life,  and  median 
life  with  the  associated  90  percent  confidence  intervals. 

Upon  completion  of  the  calculations,  the  program  calls  the 
K-S  test  and  Chi-Squared  test  subroutines  in  turn  to  provide  a  measure 
of  the  fit  of  the  Weibull  distribution  with  the  estimated  parameters 
to  the  data.  The  K-S  subroutine  "DTEST"  applies  the  Kolmogorov-Smirnov 
goodness-of-fit  test  and  prints  differences,  D,  for  each  failure  time 
in  ascending  order.  Analysis  of  the  K-S  test  is  done  by  cot.. paring  the 
D  value,  largest  in  absolute  value,  with  the  critical  value  from  an 
ordinary  table  of  K-S  D  values. 

The  Chi-Squared  test  subroutine  began  with  the  application  of 
the  Chi-Squared  test  using  a  number  of  cells  determined  by  Struges* 
rule,  K  =  1.0  +  3.3  log10  n.  However,  analysis  of  the  results  showed 
that  the  expected  and  observed  frequency  of  occurrences  were  too  small 
in  the  tails  of  the  distribution  invalidating  the  test.  The  first 
modification  to  the  subroutine  involved  combining  cells  at  the  tails 
until  the  observed  value  was  at  least  five.  This  action,  while  being 
commonly  recognized  as  being  valid,  resulted  in  a  second  problem. 
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The  Weibull  is  known  to  be  an  extreme  value,  skewed  distribution. 
Therefore,  there  is  a  small  percentage  of  area  contained  in  the  long 
right  side  tail.  The  use  of  Sturges*  rule  for  data  with  a  sample  size 
of  35  results  in  six  cells  of  equal  width.  Consequently,  when  the  data 
in  adjacent  end  cells  are  combined  to  provide  either  observed  or  expected 
frequencies  equal  to  or  greater  than  5,  the  number  of  filled  cells  may  be 
reduced  to  as  few  as  four.  Since  the  distribution  being  tested  is  the 
three-parameter  Weibull,  the  degrees  of  freedom  (k-r-1)  requires  that  the 
number  of  cells,  k,  be  at  least  five  in  order  to  have  one  degree  of 
freedom.  This  Chi-Squared  test  was  applied  to  ten  samples  of  cycles  to 
fai lure  data.  It  was  found  that  six  of  the  ten  tests  resulted  in  only 
four  filled  cells  which  provided  zero  degrees  of  freedom,  and  the  Chi- 
Squared  test  was  considered  invalid.  Thus,  it  appears  that  sample  sizes 
should  be  increased  if  the  standard  Chi -Squared  goodness-of-fit  test  is 
to  be  used. 

Another  modification  of  the  subroutine  was  made  and  tested  to  see  if 
the  sample  size  problem  could  be  circumvented.  This  modification  involved 
the  use  of  variable  cell  widths.  The  technique  described  by  Hahn  and 
Shapiro  [  4  ,  pp.  302-308]  called  for  the  calculation  of  cell  widths  to 
provide  equal  values  for  the  expected  number  of  observations.  It  was 
proposed  that  variable  cell  widths  providing  for  equal  numbers  of  actual 
observations  would  be  equally  valid  as  long  as  the  expected  frequency  was 
not  significantly  less  than  5.  A  modification  to  the  subroutine  was  made 
dividing  the  range  of  the  data  so  that  each  cell  contains  exactly  5 
failures.  For  our  sample  size  of  35  this  provides  7  cells. 

This  subroutine  was  run  for  the  same  10  sets  of  data,  and  the 
following  observations  were  made;  The  expected  frequency  for  the  seventh 
cell  was  always  less  than  5.0,  which  according  to  accepted  practices  in¬ 
validated  the  test.  It  was  observed,  however,  that  as  long  as  the  expected 
frequency  was  equal  to  or  greater  than  2.0,  the  Chi-Squared  value  did  not 
blow  up.  This  observation  was  confirmed  by  a  Monte  Carlo  simulation  of 
1,000  runs  from  which  it  was  concluded  that  Chi-Squared  errors  resulting 
from  expected  frequencies  between  2,0  and  5.0  are  insignificant. 
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A  final  modification  was  made  to  the  subroutine  using  variable 
cell  widths  and  combining  adjoining  end  cells  to  insure  that  the  expected 
number  of  observations  per  cell  equals  or  exceeds  5.  When  the  same  cycles 
to  failure  data  was  rerun  to  app.y  this  Chi-Squared  test  subroutine,  all 
tests  resulted  in  six  usable  cells  thus  providing  two  degrees  of  freedom. 

The  operation  and  accuracy  of  the  computer  program  were  verified  by 
using  the  same  input  data  used  by  Lochner  (9  )  in  his  solution  using 

Weibull  paper.  Identical  estimates  were  obtained  for  each  parameter  to 
the  degree  of  accuracy  obtainable  from  paper  plots.  The  computer  program 
provided  parameter  estimates  of  5  decimals  and  used  these  decimals  in 
subsequent  calculations.  Accuracy  of  the  D-test  subroutine  for  the 
Kolmogorov-Smirnov  goodness-of-fit  test  and  the  Chi-Squared  goodness-of- 
fit  test  were  confirmed  by  manual  (desk  calculator)  computation  of  a 
number  of  data  points.  After  validation  of  the  program  was  completed,  it 
was  used  to  reduce  the  current  data. 

3. 1 .3  SUBROUTINE  GRAPH 

The  data  reduction  programs  CYTOFR  and  WEIBULL  have  been 
extended  to  include  the  subroutine  GRAPH.  This  subroutine  performs  the 
computations  which  enable  the  Calcomp  plotter  to  draw  a  histrogram  of  the 
data  and  to  overscribe  a  curve  of  the  distribution.  The  distribution 
parameters  are  those  calculated  by  the  basic  program,  and  the  histogram 
cell  width  and  number  of  observations  are  determined  by  the  Chi-Squared  test. 
A  comparison  of  the  distribution  curve  with  the  corresponding  histogram 
provides  a  visual  measure  of  whether  the  normal,  lognormal  or  Weibull 
distribution  is  the  best  model  for  the  data.  Copies  of  subroutine  GRAPH 
are  included  in  the  program  statements  for  CYTOFR  and  WEIBULL  in 
Appendicies  C  and  D. 
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3.2  WIRE  RESEARCH  MACHINES 

3.2.1  TENSILE  DISTRIBUTIONAL  PROPERTIES 

The  results  of  the  static  strength  data  on  wires  given 
in  Section  II  are  summarized  in  Table  3.2-1.  The  static  strength  data 
was  reduced  with  the  aid  of  the  CDC-6400  program  CYTOFR.  The  normal 
distribution  was  used  to  determine  the  mean  and  standard  deviation 
parameters  for  the  ultimate  strength,  breaking  strength,  and  percent 
elongation . 

3.2.2  CYCLES-TQ-FAILURE  DISTRIBUTIONAL  PROPERTIES 

The  reduced  results  of  the  cycles-to-failure  data  on 
wires  given  in  Section  II  are  presented  in  Tables  3.2-2,  3.2-3,  3.2-4 
and  3.2-5,  and  Figures  5.2-i,,  3.2-2,  3.2-3  and  3.2-4,  for  AISI  1018, 
1038,  4130  and  4340  steels. 

The  Chi-Squazed  and  Kolmogorov-Smirnov  goodness-of-fit 
tests  were  used  to  determine  the  best-fit  distribution.  Results  of 
the  K-S  and  the  Chi-Squared  tests  are  listed  in  Table  3.2-6,  for  both 
wires  and  rods.  The  results  show  that  at  the  0.05  level  of  signifi¬ 
cance,  the  K-S  test,  does  not  reject  the  Weibull  distribution  in  22 
tests.  The  Chi-Squared  test  rejects  the  Weibull  distribution  in  7 
out  of  the  22  tests. 

The  conclusions  drawn  from  the  above  results  are: 

1.  Based  on  the  K-S  test,  the  Chi-Squared  test,  the  coefficient 
of  skewness,  ar.d  the  coefficient  of  kurtosis  values,  the 
lognormal  distribution  is  considered  generally  acceptable 
for  cycles-to-failure  data. 

2.  Based  solely  on  the  K-S  test  the  Weibull  distribution  can 
be  assumed  to  fit  cycles-to-failure  data  and  provides  an 
alternate  distribution  to  the  log  normal  (base  e)  for 
analyzing  cycles-to-failure  data. 
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3.  Based  on  the  final  Chi-Squared  test  results,  7  rejections 
out  of-22  tests,  the  Weibull  distribution  cannot  bo  consid¬ 
ered  generally  acceptable  for  cycles-to-failure  data  con¬ 
sisting  of  only  35  data  points.  Therefore,  care  must  be 
excerised  in  its  application. 


u -  «m£if 


Table  3.2-1  The  University  of  Arizona  Experimentally  Determined  Distributional  Static  Strength  Data  for  Wire. 


elongation 


Table  3.2-2  The  University  of  Arizona  Experimentally  Determined  Cycles-to-Failure  Parameters  for 

AISI  4340  SteelWire  Under  Pure  Reversed  Bending  Fatigue  Testing  Conditions  for  Design  use. 


*  Rounded  to  nearest  1,000  cycles 

*  Rounded  to  nearest  100  cycles 

'■  Beta  =  shape  parameter  or  slope 
Gamma  =  minimum  life  (cycles) 

F.ta  =  scale  parameter  (cycles) 


Table  3.2-4  The  University  of  Arizona  Experimentally  Determined  Cycles-to-Failurc  Parameters  for 

AISI  1038  Steel  Wire  Under  Pure  Reversed  Bending  Fatigue  Testing  Conditions  for  Design  use. 


Rounded  to  nearest  1,000  cycles 
Rounded  to  nearest  100  cycles 
Beta  =  shape  parameter  or  slope 
Gamma  =  minimum  life  (cycles) 

Eta  =  scale  parameter  (cycles) 


Table  3.2-5  The  University  of  Arizona  Experimentally  Determined  Cycles-to-Failure  Parameters  for 
AISI  1018  SteelWire  Under  Pure  Reversed  Bending  Fatigue  Testing  Conditions  for  Design 


0  4.0  4.5  5.0  5.5  6.0  6.5 

Log,0  CYCLES  TO  FAILURE 

Pig-  3.2-1  Statistical  S-N  Surface  for  35  specimen;  D  =  0.040  in.,  AIST  4340  steel  wire  for  results 
summarized  in  Tables  3.2-1  and  3.2-2. 


0  4.0  4.5  5.0  5.5  6.0  6.5  7. 

LogIO  CYCLES  TO  FAILURE 

Fig.  3.2-2  Statistical  S-N  surface  for  35  specimens;  D  =  0.040  in.,  AISI  4130  steel  wire  for  results 
summarized  in  Tables  3.2-1  and  3.2-3. 


Loglo  CYCLES  TO  FAILURE 

Fig.  3.2-3  Statistical  S-N  Surface  for  35  specimen;  D  =  0.040  in.,  AISI  1038  steel  wire  for  results 
summarized  in  Tables  3.2-1  and  3.2-4. 


0  4.0  4.5  5.0  5.5  6.0  6.5  1 

Log,0  CYCLES  TO  FAILURE 

Fig.  3.2-4  Statistical  S-N  Surface  for  35  Specimen;  D  =  0,040  in.,  AISI  1018  steel  wire  for  results 
summarized  in  Tables  3.2-1  and  3.2-5. 


Table  3.2-6  K-S  and  Chi-Squared  test  results  for  Weibuil  distribution 


Material  Stress  Level 

AISI  - 

psi 

Steel 

Wire  D  =■  0.G40  Stress  Ratio  =  » 


4130 

4130 

4130 

4130 

4130 

1038 

1038 

1038 

1038 

1018 

1018 

1018 

4340 

4340 

4340 

4340 

4340 


0.125 

0.375 


67.700 
70,000 

72.500 

74.700 

77.800 

64.500 

67.200 

69.200 
72,300 

57.200 
60,000 

62.800 

73.500 

78.100 

80.100 

84.700 
90,000 


d  =  0.0937 
Stress  Ratio 


Maximum  D 
Value  (K-S) 

dof 

0.078 

2 

0.124 

2 

0.066 

2 

.  0.114 

2 

0.112 

3 

0.102 

2 

0.141 

2 

0.167 

3 

0.097 

3 

0.114 

2 

0.107 

2 

0.064 

1 

0.134 

3 

0.106 

3 

0.328 

3 

0.059 

3 

0.102 

2 

rr  oa 

.106 

3 

.107 

2 

.120 

2 

0.250  d  =  0.0937 

0.375  Stress  Ratio  =  " 


Critical  Value  at  .05  Level  of  Significance: 
K-S  D  Value  0.224. 

Chi-Squared,  1  dof  3,841. 

Chi-Squared,  2  dof  5.991. 

Chi-Squared,  2  dof  7.815. 


Chi -Squared 
Value 


2.127 

8.102 

3.674 

4.097 

3.390 

3.708 

6.507 

4.245 

3.863 

3.724 

4.116 

2.546 

3.075 

4.993 

69.373 

1.227 

7.581 
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3.2.3  STAIRCASE  TEST  AND  DISTRIBUTIONAL  PROPERTIES 

The  staircase  fatigue  tests  of  0.040  in.  dia.  wire  for  steels  of 
AISI  4130,  4340,  1038  and  1018  are  being  run  and  shall  be  reported  in 
ONR  VI  with  completed  S-N  diagrams. 

Knowing  that  the  objective  of  staircase  fatigue  testing  is  to 
determine  the  mean  and  the  standard  deviation  of  the  stress  to  failure 
for  various  metals  with  various  geometries,  a  standard  procedure  for 
obtaining  this  information  was  established  and  appears  next. 

HOW  MANY  SPECIMENS  TO  RUN? 

In  determining  the  number  of  test  specimens  to  use  one  must  account 
for  the  (1)  cost,  (2)  time  and.  (3)  the  sufficiency  of  the  sample  size  to 
draw  significant  conclusions  from  the  results. 

Since,  staircase  testing  is  concentrated  near  the  mean  fatigue 
strength  the  number  of  specimens  tested  may  be  less  than  for  the  "Probit" 
method  Qo,  p.  12],  which  gives  results  for  a  wider  range  of  stress  values. 

In  general  it  is  recommended  that  at  least  30  specimens  be  tested  fio,  p.13]. 
Thirty  specimens  also  correspond  to  the  minimum  number  of  specimens  needed 
for  determining  95  percent  confidence  intervals  for  a  population  standard 
deviation,  a.  For  the  data  collected  through  staircase  testing  to  be 
meaningful  it  is  recommended  that  at  least  35  specimens  be  tested. 

HOW  10  SELECT  THE  NUMBER  OF  CYCLES  A  SPECIMEN'  SHOULD  RUN. 

The  number  of  cycles  to  run  is  determined  by  the  cycles  of  life  for 
which  the  stress-to-failure  distribution  is  needed.  For  the  endurance 
stress-to-failure  (strength)  distribution  it  must  be  insured  that  the 
specimen  life  chosen  is  beyond  the  knee  of  the  S-N  diagram.  It  has  been 
found  that  normally  the  knee  of  the  S-N  curve  will  occur  between  10  and 

7 

10  cycles  for  steel  specimens.  It  is  recommended  here  that  the  specimens 
be  tested  for  3  x  10  cycles  or  to  failure  whichever  occurs  first.  Once 
the  specimen  has  reached  3  x  10^  cycles  without  failure  the  testing  of  that 
particular  specimen  should  be  terminated,  and  a  fresh  specimen  be  tested 
at  one  stress  increment  higher. 


v 


.a— i 
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HOW  TO  DETERMINE  THE  INITIAL  STRESS  LEVEL  FOR  UNGROOVED  SPECIMENS. 


When  initiating  a  staircase  fatigue  test  it  is  desirable  to  begin 
the  test  at  or  near  the  mean  endurance  strength  of  the  material  under 
investigation.  It  has  been  determined  experimentally  that  a  good  estimate 
for  the  beginning  stress  level  of  a  staircase  test  is  0.54  of  the  material's 
ultimate  strength.  Therefore,  when  initiating  a  fatigue  staircase  test  of 
ungrooved  specimens  the  following  initial  stress  level  is  recommrnded:  _ 


S  =  0.54  Sm4. .  . 

st  Ultimate 


(3.2-1) 


HOW  TO  DETERMINE  THE  INITIAL  STRESS  LEVEL  FOR  GROOVED  SPECIMENS. 


When  fatigue  testing  grooved  steel  specimens  the  fatigue  stress 
concentration  due  to  the  presence  of  the  groove  must  be  taken  into  account. 
To  do  this  a  factor,  X^,  known  as  the  fatigue  stress  concentration  factor 
must  be  taken  into  account.  A  summary  of  estimated  fatigue  stress  con¬ 
centration  factors  is  given  in  Chapter  2.2.1.  Once  the 

value  for  the  groove  configuration  has  been  determined  the  following 
relationship  is  recommended  for  the  starting  stress  of  a  staircase  test 
involving  grooved  specimens: 


0.54  S, 


Ultimate 


(3.2-2) 


HOW  TO  DETERMINE  STRESS  STEPS. 


When  conducting  the  staircase  endurance  tests  on  the  modified 
Wiedemann  and  the  unmodified  Wiedemann  research  machines  it  is  suggested 
that  stress  increments  of  2,000  psi  be  used.  The  reason  for  this  is  that 
stress  levels  of  less  than  2,000  psi  are  extremely  difficult  to  obtain  due 
to  the  fact  that  only  about  0.05  pounds  of  pan  weight  is  required  to  change 
the  stress  level  by  1,000  psi.  Therefore,  it  is  recommended  that  stress 
increments  of  2,000  psi  be  used. 

For  the  wire  fatigue  research  machines  the  stress  increment  is 
determined  from  the  calibration  curve  by  selecting  the  change  in  stress 
level  corresponding  to  0.5°  change  in  bend  angle.  While  this  procedure 
results  in  different  stress  increments  for  different  test  groups,  it  provides 
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v—  for  greater  consistency  in  setting  the  bend  angle  carriage.  The  calibration 

curve  is  quite  linear  over  the  endurance  range;  therefore,  use  of  a  0.5 
degree  change  in  bend  angle  results  in  linear  stress  increments. 


HOW  TO  DETERMINE  THE  STRESS  LEVEL  TO  USE  FOR  THE  SECOND  AND 
SUBSEQUENT  SPECIMEN'S. 

g 

If  the  first  test  specimen  fails  before  3  x  10  cycles  then  a  second 
test  specimen  should  be  tested  at  a  stress  level  2,000  psi  below  the  first 
stress  level.  If  the  first  test  specimen  survives  to  3  x  10^  cycles  the 
second  specimen  should  be  tested  at  a  stress  level  2,000  psi  higher  than 
the  first  stress  level,  as  shown  in  Fig.  3.2-5  •  This  procedure  should 

be  followed  for  all  subsequent  specimens  tested. 

WHEN  TO  TERMINATE  A  STAIRCASE  TEST  AND  WHICH  DATA  TO  USE. 

In  staircase  testing  the  first  few  specimens  are  not  included  in 
the  test  if  no  change  of  mode  (from  failure  to  success  or  vice  versa)  occurs. 
The  test  is  considered  to  have  started  with  the  specimen  prior  to  the  one 
where  the  first  change  of  mode  occurs.  Only  about  half  of  the  test  results 
are  used  in  computing  the  mean  and  standard  deviation  of  the  strength  dis¬ 
tribution.  Either  the  successes  or  the  failures  are  used,  whichever  are 
the  fewer.  Therefore,  33  specimens  should  be  tested  after  the  first  mode 
change  occurs.  Thus,  the  staircase  test  may  be  terminated  when  33  specimens 
have  been  tested  after  the  initial  mode  change  has  occurred  for  a  total  of 
35  good  specimens. 

HOW  TO  PLOT  THE  STAIRCASE  DATA. 

Standard  forms  have  been  developed  for  plotting  staircase  test  data. 
On  the  ordinate  of  the  form  should  appear  the  nominal  stress  level  in  Kpsi, 
and  on  the  abscissa  the  number  of  specimens  tested.  To  plot  the  data  a 
legend  of  black  dots  and  white  dots  are  used.  A  black  dot  signifies  that 
the  test  specimen  failed  prior  to  3  x  10^  cycles  and  a  white  dot  signifies 
that  the  specimen  survived  bo  3  x  10^  cycles.  All  data  should  then  be 
plotted  accordingly,  as  shown  in  Fig.  3.2-5.  • 


HOW  TO  CALCULATE  THE  STRESS -TO-FAILURE  DISTRIBUTIONS  MEAN  AND 
STANTARD  DEVIATION  ASSUMING  A  NORMAL  DISTRIBUTION. 


106 


The  equations  given  here  assume  that  the  strength  distribution  is 
normal.  Some  distributions  can  be  made  normal  by  simple  transforms.,  such 
as  the  lognormal.  If  the  data  can  be  assumed  normal  or  transformed  to 
normal,  the  statistical  analysis  of  the  data  is  as  follows  (n  ,  p.  114]: 

The  estimate  for  the  mean  strength  is  given  by 

X  =  Xq  +  d  l,VN  +  1/2],  (3.2-5) 

where  X  =  the  estimate  of  mean  strength,  psi. 

Xq  =  the  lowest  stress  level  used,  psi. 

(If  the  fewer  events  are  the  successes  then  the  lowest 
stress  level  used  is  the  lowest  level  of  the  staircase 
plot.  In  Fig.  3.2-5  this  level  is  the  23.000  psi  level. 

If  the  fewer  events  are  the  failures  then  the  lowest  stress 
level  used  is  the  lowest  level  in  the  staircase  plot  where 
failures  have  occurred.  In  Fig.  3.2-5  this  level  is  the 
25,000  psi  level. 

d  =  the  stress  increment,  psi. 

(Normally  use  2,000  psi.) 

N  =  total  number  of  successes  or  failures  whichever  are  the 
fewer  in  number. 

The  quantity,  A,  in  Eq.(5.2-3.)  is  defined  as 

H 

A  =  il0  inr  (3.2-4) 

In  Eq, (3.2-4)  "i"  is  the  numerical  order  number  of  each  stress  level 
with  the  lowest  stress  level  used  in  the  calculations  having  i  =  o,  the 
next  stress  level  having  i  =  1,  etc.,  and  n^  is  the  number  of  failures  or 
successes,  whichever  has  been  chosen  as  the  basis  for  the  calculations, 
occurring  at  each  of  the  "i"  stress  levels.  For  an  example  of  the 
calculations  see  Table  3.2-7  .  H  is  the  order  number  of  the  highest 
stress  level  used. 


lx  |x 
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Table  3.2-7  Reduced  Data  for  Group  163,  AISI  1038  Steel.  Wiedemann 
Fatigue  Machine. 

6 

Staircase  method  at  3  x  10  cycles 

Number  of  Useful  Specimens:  35 

Specimen  geometry:  D  =  0.375  in. 

d  =  0.2700  in. 
r  =  0.062  in. 


Alternating 

Stress 

psi 

27,000 

25,000 

23,000 


n. 

1 

Successes 

in. 

X 

.2 

1  n. 

1 

9 

18 

36 

5 

5 

5 

3 

0 

0 

N  *  17 

A  =  23 

B  =  41 

d  =  stress  increment  =  2,000  psi 
Xq  =  lowest  stress  level  =  23,000  psi 
=  mean  (estimate) 

=  Xq  +  d  A/N  +  1/2  =  23,000  +  2,000 
X  =  26,706  psi 

s  =  standard  deviation  (estimate) 
s  =  1.620  d  (NB-A2)/N2  +0.029  =  1.620(2 

s  =  1,976  psi 


,000)  Ifllsiii *  0.029] 

17 


- 2.  S^JL - 
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In  Eq.  (3.2-3)  the  positive  (+)  sign  before  the  one-half  (1/2) 
is  used  if  the  less  frequer.t  events  are  the  successes,  and  the  negative 
(-)  sign  if  the  failures  are  the  less  frequent. 

The  estimate  of  the  standard  deviation  of  the  strength  distri¬ 
bution  is  given  by 

.«  =  1.620  d  [(NB  -  A2)/N2  +  0.029],  (3.2-5) 

where  S  =  standard  deviation  estimate,  psi. 

The  dimensionless  quantity  B  is  defined  as  follows: 

H  2 

B  =  Z  in..  (3.2-6) 

i 

1=0 

For  an  example  of  how  to  apply  Eqs.  (3.2-5)  and  (3.2-6)  see  Table  3.2-7 
and  Fig.  3.2-5. 

See  Appendix  E  for  a  PDP-8  Focal  computer  program  to  calcui ate 
x  and  S  for  the  staircase  test. 


,  -l  x.  -i.  .fc^.  ^ 
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RECOteiENDATIONS 

1.  It  is  recommended  that  at  least  35  specimens  he  tested  when  conducting 
staircase  tests. 

2.  It  is  recommended  that  3  x  10^  cycles  be  used  as  the  cut-off  point 
when  conducting  staircase  endurance  tests. 

3.  When  determining  the  starting  stress  level  for  ungrooved  specimens  use 
the  following  relationship: 


=  0.54  S 


Ultimate' 


(3.2-1) 


4. 


When  determing  the  starting  stress  level  for  grooved  specimens  use 
the  following  relationship: 


0.54  S 

Ultimate 


(3.2-2) 


5.  Stress  increments  of  2,000  psi  are  recommended  when  conducting 
staircase  tests. 


6.  It  is  recommended  that  a  staircase  test  not  be  terminated  until  33 
specimens  have  been  tested  after  the  first  mode  change  has  occurred, 
for  a  total  of  35  useful  specimens. 


7.  If  when  reducing  the  data  the  failures  are  being  used  the  lowest  stress 
level  should  be  excluded  from  the  calculations. 

8.  If  when  reducing  the  data  the  successes  are  being  used  the  highest  stress 
level  should  be  excluded  from  the  calculations. 

9.  When  reducing  the  data  the  basis  of  the  calculations  should  be  the  less 
frequent  event;  i.e.,  if  there  are  fewer  failures  than  successes  then 
the  failures  should  be  used,  if  there  are  fewer  successes  than  failures 
then  the  successes  should  be  used. 

10.  It  is  recommended  that  the  following  relationship  be  used  to  calculate 
the  mean  of  the  strength  distribution: 
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X  =  Xo  +  d  (A/N  +  1/2)  .  (3.2-3) 

11.  It  is  recommended  that  the  following  relationship  be  used  to 
calculate  the  standard  deviation  of  the  strength  distribution: 

S  =  1.620  d  [(NB  -  A2)  /  N2  + 


0.029]  . 


(3.2-5) 
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3.3  WIEDEMANN  RESEARCH  MACHINES 

3.5.1  TENSILE  DISTRIBUTIONAL  PROPERTIES 

Table  3.3-1  presents  the  experimental  data  generated  to 
determine  the  distributional  static  strength  data  for  ungrooved  steel 
rod  under  the  ONR  Contract.  Listed  in  the  table  are  the  type  of 
material  tested,  the  geometry,  the  composition,  sample  size,  yield 
strength,  ultimate  strength,  breaking  strength,  percent  elongation  and 
their  parameters. 

3.3.2  CYCLES-T0-FA1LURE  DISTRIBUTIONAL  PROPERTIES 

Figures  3.3-1  and  3.3-2  present  two  S-N  diagrams 
generated  for  AISI  4130  grooved  steel  rod  in  reversed  bending  for  r  = 
0.250"  and  r  =  0.125",  respectively.  The  experimental  data  for  these 
figures  are  given  in  Tables  3.3-2  and  3.3-3,  respectively.  Listed  in 
each  table  are  the  type  of  material  tested,  the  geometry,  and  the  normal 
and  lognormal  distribution  parameters  of  the  cycles-to-failure  data  at 
various  alternating  stress  levels.  These  tables  provide  a  quick  means 
for  obtaining  the  mean  and  standard  deviation  of  the  cycles-to-failure 
distribution  at  various  alternating  stress  levels.  The  S-N  diagrams 
appearing  in  Figs.  3.3-1  and  3.3-2  contain  the  mean  line  of  the  cycles- 
to-failure  distributions  and  the  ±  3  o  envelopes  using  the  lognormal 
distribution  as  the  one  best  representing  the  cycles-to-failure  data 
at  various  stress  levels. 


3.3.3  STAIRCASE  TEST  AND  DISTRIBUTIONAL  PROPERTIES 

Table  3.3-4  presents  the  experimental  data  generated  to 
determine  the  distributional  parameters  for  staircase  testing  of 
AISI  1038  steel  rod  under  pure  reversed  bending.  Listed  in  the  table 
are  the  group  number,  geometry,  number  of  useful  specimens,  and  the 
mean  and  standard  deviation  of  the  endurance  strength  distribution  for 
3  x  10^  cycles  of  life. 
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Table  3.3-1  The  University  of  Arizona  experimentally  determined  distributional 
static  strength  data  for  ungrooved  steel  rod. 


AISI  4130  Steel 

AISI  4340  Steel 

AISI  1038  Steel 

Composition 

.S05C  . S3Mn  .012P 

.385C  . 72Mn  .009P 

.39C  . 80Mn 

.OllSu  . 22Mo  .2751 

.OllSu  .28Si  .78Cr 

.018Su  .013P 

.98Cr  . 15Ni  .16Cu 

. 175Ni  . 26Mo 

Diameter 

in. 

.250 

.270 

.250 

Yield 

Mean* 

87 

Strength 

Std.Dev.** 

2.0 

Kpsi 

C.  V.*** 

.023 

Ultimate 

Mean 

106 

116 

69 

Strength 

Std.Dev. 

2.3 

1.3 

1.3 

Kpsi 

C.  V. 

.020 

.011 

.019 

Breaking 

Mean 

152 

158 

Strength 

Std.Dev. 

4.9 

6.4 

Kpsi 

C.  V. 

.032 

.040 

Percent 

Mean 

8.5 

20.9 

Elongation 

Std.Dev. 

0.9 

2.3 

c.v.**** 

0.10 

0.11 

Sample 

Size 

35 

35 

35 

Source 

4,  pp.  471 

,473,475 

2.  p>  205 

4,  p.  488 

*  Rounded  to  nearest  1  Kpsi 

**  Rounded  to  nearest  0.1  Kpsi 

***  Coefficient  of  variation  =  — 

S 

ani 

****  Coefficient  of  variation  =  — — 

EL 


EL 


Elongation 
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0  3.0  3.5  4.0  4.5  5.0  5.5  6.0  * 

Log,0  CYCLES  TO  FAILURE 

Fig.  3.3-2  Experimentally  determined  S-N  diagram  using  ONR  research  results  for  AISI  4330  grooved  steel 
rod  in  reversed  bending,  D  =  0.375",  d  =  0.0937"  and  r  =  0.125". 


Table  3.3-2  The  University  of  Arizona  Experimentally  Determined  Cycles -to-Fai lure  Parameters  for 

AISI  4130  Steel  Rod  Under  Pure  Reversed  Bending  Fatigue  Testing  Conditions  for  Design  use. 


Rounded  to  nearest  1,000  cycles 
Rounded  to  nearest  100  cycles 
Beta  =  shape  parameter  or  slope 
Gamma  =  minimum  life  (cycles) 

Eta  =  scale  parameter  (cycles) 


Table  3.3-3  The  University  of  Arizona  Experimentally  Determined  Cycles-to-Failure  Parameters  for 

AISI  4130  Steel  Rod  Under  Pure  Reversed  Bending  Fatigue  Testing  Conditions  for  Design  use. 
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Table  5.3-4 


The  University  of  Arizona  experimentally  determined 
staircase  parameters  for  AISI  1058  steel  rod  under  Dure 
reversed  bending,  fatigue  testing  conditions  for  3  x  10b 
cycles  of  life. 


Group 

No/ 


Geometry 

in. 


r  -  .031 
d  =  .2700 
D  =  .375 


r  =  .062 
d  =  .2700 
D  =  .375 


.125 

.2700 

.375 


Number  of 
Useful  Specimens 


Stress  (psi) 
Mean*  Std.  Dev.* 


.250 

.2700 

.375 


1.87 

.2700 

.375 


33 

23,800 

400 

36 

26,700 

2,000 

37 

30.,  444 

1,500 

35 

32,200 

1,200 

36 

j 

i 

37,100 

2,100 

Rounded  off  to  nearest  100  psi . 


j 


118 

3.3.4  METHODOLOGY  FOR  CONSTRUCTING  S-N  DIAGRAMS  EMPIRICALLY  FROM 

STATIC  ULTIMATE  AND  ENDURANCE  STRENGTH  DATA  FOR  GROOVED  AND 

UNGROOVED  SPECIMENS. 

Gne  of  the  objectives  of  fatigue  research  is  to  obtain 
expressions  to  estimate  the  life  of  a  machine  part,  manufactured  from  a 
given  material  and  subjected  to  cyclic  loading. 

The  fatigue  strength  of  a  material  may  be  defined  as  the 
maximum  stress  that  can  be  applied  repeatedly  to  the  material  without 
causing  failure  in  less  than  a  certain  finite  number  of  cycles.  The 
endurance  strength  is  the  maximum  stress  that  can  be  applied  repeatedly 
to  a  material  for  a  large  number  of  cycles  (beyond  the  knee  of  the  S-N 
diagram  or  for  10^  cycles  or  more  for  steels)  without  causing  failure. 

The  relationship  between  the  magnitude  of  repeatedly  applied  stress  and 
the  number  of  cycles  to  failure  is  conventionally  presented  in  the  graphical 
form  known  as  the  S-N  diagram  [l2,PP*  1-2]. 

Empirical  results  have  shown  that  fatigue  life  is  an  exponential 
function  of  the  maximum  level  of  repeated  stress.  Other  empirical  relations 
developed  for  strain-cycle  fatigue  show  that  cyclic  life  is  also  related  to 
the  applied  strain  range  exponentially.  Therefore,  S-N  diagrams  are 
usually  presented  in  semi-log  or  log-log  form  [l2,p.2];  however,  the  log- 
log  form  is  the  most  prefered. 

As  the  fatigue  properties  of  a  material  are  determined  by  a 
variety  of  factors,  the  shape  of  the  S-N  curve  varies  according  to  the 
conditions  represented.  Factors  influencing  the  shape  of  the  S-N  curve 
are  the  type  of  loading,  temperature,  operating  environment,  surface  finish, 
stress  concentration  and  size  effects.  In  this  discussion  machine  parts 
subjected  to  fully-reversed  bending  stresses  are  investigated. 

TYPE  OF  DESIGN  DATA  NEEDED  FOR  CONSTRUCTING  THE  S-N  DIAGRAM 

There  are  a  number  of  suggested  methods  for  constructing  S-N  diagrams, 
the  most  accurate  being  that  based  on  experimental  results.  However,  this 
is  also  the  most  costly  and  time  consuming.  The  ultimate  goal  here  is  to 
develop  a  method  for  empirically  constructing  acceptable  S-N  diagrams  from 
minimum  and  easily  obtainable  data,  thereby  saving  both  time  and  money. 


•  1  -|  i-  - - -  Mlttfc 
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The  experimentally  determined  S-N  diagram  is  based  upon  much  testing. 

The  static  ultimate  tensile  strength  for  both  grooved  and  ungrooved 
specimens  must  be  established,  as  well  as  the  endurance  strength.  Then 
at  intermediate  stress  levels  between  the  static  ultimate  and  the  endurance 
strength,  the  cycles-to-failure  or  stress-to-£ailure  distributions  must 
be  generated.  This  requires  over  200  specimens  and  extensive  testing  time. 
Therefore,  it  is  highly  desirable  to  establish  empirical  relationships  to 
give  good  estimates  of  the  S-N  diagrams  based  on  a  minimun  amount  of  data. 

Most  available  empirical  methods  for  constructing  the  S-N  diagram  use 
the  ungrooved  static  ultimate  tensile  strength  and  make  certain  modifications 
to  estimate  the  endurance  strength.  This  requires  the  determination  of  the 
static  ultimate  tensile  strength  and  certain  modifications  to  estimate 
the  endurance  strength.  This  requires  the  determination  of  the  static 
ultimate  tensile  strength  of  ungrooved  specimens  only,  and  modifying  factors. 
However,  there  is  one  major  pitfall  involved  in  this.  Preliminary  test 
data  generated  by  ONR  research,  indicates  that  the  static  ultimate  tensile 
strength  of  grooved  specimens  is  substantially  greater  than  the  static 
ultimate  tensile  strength  of  ungrooved  specimens.  Thus  using  ultimate 
strength  data  of  ungrooved  specimens  to  construct  an  S-N  diagram  for  grooved 
specimens  tends  to  yield  very  conservative  results  at  the  low-cycle  end  of 
of  the  S-N  diagram.  A  full  test  plan  for  investigating  the  effects  of  using 
grooved  ultimate  tensile  strength  in  constructing  S-N  diagrams  is  now  being 
planned. 


EMPIRICAL  PROCEDURE  FOR  CONSTRUCTING  S-N  DIAGRAMS  FOR  UNGROOVED 
SPECIMENS 

The  method  for  constructing  the  S-N  diagram  for  ungrooved  specimens 
is  considerably  simpler  than  for  grooved  specimens.  Shigley  [6»  pp. 

160-191]  and  Juvinall  [7,  pp.  210-224]  recommend  identical  methods.  In 
the  absence  of  specific  endurance  strength  data  they  both  recommend  that 
the  endurance  strength  be  estimated  as  [6,  p.162],  [7,  p.  211] 


S  =  0.5  S  , 

e  u 


(3.3-1) 


i  I'r.'.ttf 
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;i 

i-i 

where  I 


S  =  estimated  endurance  strength,  psi.  I 

e  ,j 

Su  =  ungrooved  static  ultimate  tensile  strength,  psi. 

i 

To  estimate  the  stress  level  at  10  cycles  on  the  S-N  diagrams  both  4 

Shigley  [6,  p.  163]  and  Juvinall  [7,  p.  211]  recommend  the  following  rule 
of  thumb:  .J 


S,n3  .  =  0.9  S  , 

10  cycles  u 


(3.5-2) 

i 


where 

3 

S, „3  ,  =  alternating  stress  at  10  cycles,  psi. 

10  cycles 

,  i 

PUBLISHED  METHODOLOGIES  FOR  CONSTRUCTING  S-N  DIAGRAMS  FOR 

GROOVED  SPECIMENS  j 

One  published  method  for  constructing  S-N  diagrams  is  that  recommended 
by  Shigley  [6,  pp.  161-199].  Shigley  states  the  following  [6,  pp.  161-  ! 

162]: 

"Great  numbers  of  rotating  beam  tests  have  shorn  that  the 
endurance  limit  varies  from  40  to  60  percent  of  the 
ultimate  tensile  strength  for  steels  with  an  ultimate 
strength  up  to  200,000  psi.  It  has  become  a  more  or  less 

standard  practice,  when  the  results  of  tests  are  not  1  j 

available,  to  estimate  the  endurance  limit  for  steels  as 

j 

S  '  =  0.5  Sy  for  Su  ^  200,000  psi,  ■  1 

and 

S'  =  100,000  psi  for  Sy  >  200,000  psi."  (3.3-3)  J 


Shigley  then  suggests  modification  of  Eq.  (3.3-1)  to  consider  the 
effects  of  surface  finish,  size  and  stress  concentration^,  p.  166], 

The  modified  relationship  is 


S 

e 


=  k  k.  k 
a  b  e 


(0.5  S  1 
u  > 


(3.3-4) 


«■** 


where 
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Se  =  corrected  endurance  limit,  psi 

k  =  surface  factor 

a 

k,  =  size  factor 

b 

kg  =  modifying  factor  for  stress  concentration  (This  is 

not  the  same  as  the  stress  concentration  factor,  K  ). 

Shigley  discusses  the  effects  of  surface  finish  on  fatigue  life  and 
presents  Fig.  3.3-3  '  for  obtaining  the  surface  factor,  k  ,  for  steels 

ci 

when  the  tensile  strength  and  the  surface  finish  are  specified.  The 
surface  factor  for  wrought  and  cast  aluminum,  magnesium  and  other  non- 
ferrous  materials  may  be  taken  as  unity,  because  the  endurance  limits 
listed  by  the  manufacturers  of  these  materials  usually  include  the  effect 
of  surface  finish  [6,  pp.  166-167] . 

It  is  clear,  however,  that  the  surface  finish,  at  least  for  steel 
parts,  has  a  significant  effect  on  the  endurance  limit.  Fig.  3.3-3 
shows  that  the  endurance  limits  for  hot-rolled  and  as-forged  parts  improve 
very  little  with  increase  in  tensile  strength.  It  also  shows  that,  when 
fatigue  loading  is  an  important  design  factor,  there  is  little  to  be 
gained,  even  with  a  ground  finish,  by  selecting  steels  having  tensile 
strengths  greater  than  200,000  psi  [6>  p.  166], 

The  standard  rotating  beam  test  gives  the  endurance  limit  for  a 
specimen  of  0.30  in.  in  diameter  at  the  test  section.  When  specimens  of 
larger  diameter  are  tested  subjected  to  completely  reversed  stress  in 
bending  or  torsion,  it  is  found  that  the  endurance  strength  is  10  to  15 
percent  lower  for  specimens  up  to  2  in.  in  diameter  [6,  p.  168],  Therefore, 
Shigley  suggests  a  size  factor,  k^,  of  0.85  for  specimens  with  a  diameter 
greater  than  0.30  in.  subjected  to  reversed  bending  or  torsional  stresses. 

No  reduction  was  found  for  completely  reversed  axial  loading;  consequently 
k^  =  1.0  for  axial  loads  fe ,  p.  168],  Shigley  then  suggests  that  for 
machine  elements  larger  than  2  in.  in  diameter,  adjustments  somewhat  larger 
than  10  percent  be  made  to  the  endurance  strength  but  specific  values  are 
not  given. 
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As  a  fatigue  failure  almost  always  originates  at  a  discontinuity 
the  effect  of  stress  concentration  should  be  brought  in.  Usually  the 
crack  begins  at  a  groove,  a  notch,  a  shoulder  or  at  the  edge  of  a  hole, 
but  it  may  also  start  at  a  tool  mark,  a  scratch  or  even  at  an  inspection 
mark.  In  the  investigation  of  fatigue  failures  it  is  found  that  some 
materials  are  much  more  sensitive  to  notches  than  others.  This  means 
that  the  stresses  for  insensitive  materials  need  not  be  increased  by  the 
full  amount  of  the  geometric  stress  concentration  factor,  For  this 
reason  it  is  useful  to  define  a  fatigue  stress-concentration  factor,  K^., 
as  follows  [6,  p.  170]: 


K  _  Endurance  limit  of  ungrooved  specimes 

f  Endurance  limit  of  grooved  specimens 


(3.3-5) 


Shiglev  then  denotes  kg  as  a  fatigue-strength  reduction  factor  and  defines 
this  modifying  factor  for  fatigue  stress-concentration  as 


(3.3-6) 


In  order  to  account  for  the  various  sensitivities  of  materials  to  notches 
a  factor,  q,  known  as  notch  sensitivity  is  introduced  as 


(3.3-7) 


where  q  is  between  zero  and  unity.  Equation  (3.3-7)  shows  that  if  q  =  0, 

=  1  and  the  material  has  no  sensitivity  to  notches  at  all.  On  the  other 
hand,  if  q  =  1,  then  =  K  and  the  material  has  full  notch  sensitivity. 

In  analysis  for  design  work  one  first  determines  K  from  the  geometry  of  the 
part.  Then,  the  material  having  been  specified,  q  can  be  found.  Equation 
(3.3-7  )  also  provides  a  means  for  estimating  without  generating  test 
data.  Solving  Eq.  (3.3-7)  for  yields 


Kf  =  1  +  q  (Kt  -  1)  . 


(3.3-8) 
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Charts  prepared  by  Peterson  f5,  pp.  9-10]  give  notch  sensitivity  values 
for  steels  in  reversed  bending  which  have  been  used  widely.  Thus,  after 
locating  the  notch  sensitivity  value  and  the  theoretical  stress  concen¬ 
tration  value  from  other  charts  given  by  Peterson  [5,  pp.  47-50]  for  the 
machine  part  under  investigation,  Eq. (3.3-8  )  can  be  used  to  obtain  the 
value. 

Shigley  suggests  that  whenever  there  is  any  doubt  concerning  the 
proper  notch  sensitivity  value  to  use,  one  can  always  make  and 

err  on  the  safe  side.  Shigley  also  suggests  that  when  the  groove  radius 
is  quite  large  is  not  far  from  unity  and  the  error  of  assuming 
equal  to  K^_  will  be  quite  small  [6,  p.  172]. 

USE  OF  SHIGLEY1 S  METHOD  TO  CONSTRUCT  AN  S-N  DIAGRAM  FOR  ONR 
GROOVED  TEST  SPECIMENS 

Shigley' s  method  for  constructing  an  S-N  diagram  will  be  applied  to 
parts  like  the  test  specimen  used  in  this  research.  The  results  will  then 
be  compared  with  the  S-N  diagram  obtained  by  the  use  of  the  experimental 
data  generated  during  this  research  with  specimens  of  AISI  4130  steel  rod 
subjected  to  fully  reversed  bending  and  having  an  outside  diameter  of 
D  =  0.375  in.,  a  diameter  at  the  base  of  the  groove  of  d  =  0.0937  in., 

and  a  groove  radius  of  r  =  0.250  in.  The  ungrooved  ultimate  tensile 
strength,  Su,  was  found  experimentally  to  have  a  mean  value  of  106,000  psi. 
The  groove  in  the  specimen  was  machined  and  then  polished  to  approximately 
4^-in.  per  in.  The  polishing  technique  involved  the  use  of  medium  grit 
emery  paper  and  then  a  fine  grit  lapping  compound. 

Using  Shigley' s  approximation  for  5  1  as  given  by  Eq.  (3.3-3) 

S'  =  0.5  (106,000), 
or 

S  r  =  53,000  psi. 


The  modifying  factors  are: 
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kb  =  0.85  [6,  p.  168], 
k  =  0.89  [6,  p.  167], 

cL 

and 

v  =  =  _ i _ 

Ke  Kf  1  +  q(Kt  -  1)  * 


where  q  values  are  obtained  from  Peterson’s  notch  sensitivity  charts 


Substitution  of  106,000  psi,  the  static  ultimate  tensile  strength,  into 
Eq.  (3.3-10)  gives 


^10°  cycles 


0.9(106,000) 


or 

S, -3  .  =  95,400  psi. 

10  cycles  *  r 

Thus  the  estimated  S-N  diagram  can  now  be  constructed  for  the  specimen 
involved,  using  Shigley's  recommendations,  as  shown  in  Fig.  5.3-3. 


COMPARISON  OF  SHIGLEY'S  EMPIRICAL  METHOD  WITH  EXPERIMENTAL  RESULTS 

Upon  comparing  Fig. 3, 3.4  (S-N  diagram  constructed  using  Shigley's 

methods)  with  Fig.  3.3-1  ,  which  was  constructed  using  exclusively 

experimental  data  generated  during  this  research,  it  is  found  that  Shigley's 

3 

method  gives  very  conservative  results.  The  mean  value  at  10  cycles  by 

Shigley's  method  is  95.4  Kpsi,  whereas  in  Fig.  3.3-1  the  value  is  112  Kpsi.  It 
also  gives  conservative  results  for  the  endurance  strength  at  10  cycles. 

With  Shigley's  method  the  endurance  strength  is  38.1  Kpsi  while  the  ONR 
experimental  results  give  a  value  of  54  Kpsi.  As  can  be  seen,  Shigley's 
method  is  more  conservative  toward  the  endurance  end  of  the  S-N  diagram  and 
the  entire  curve  generated  by  Shigley's  method  falls  outside  the  -3  0  limit 
at  both  ends  of  the  experimentally  determined  S-N  diagram.  Although 
Shigley's  method  is  very  easy  to  apply,  as  it  requires  the  knowledge  of  the 
static  ultimate  strength  only,  it  should  be  employed  when  data  such  as  is 
being  generated  by  this  research  is  not  available.  Even  then  being  so  far 
on  the  conservative  side  has  its  penalties  in  the  form  of  cost  and  weight, 
and  it  ignores  the  distributional,  statistical  nature  of  the  fatigue 
phenomenon. 
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JUVINALLS 1  METHOD  FOR  CONSTRUCTING  S-N  DIAGRAMS  FOR  GROOVED 
SPECIMENS  IN  BENDING 


The  S-N  diagram  discussed  by  Juvinall  [7,  pp.  210-262]  is  constructed 
in  basically  the  same  manner  as  Shigley's  [6  .  pp.  160-191].  Both  methods 
base  the  entire  procedure  on  the  use  of  the  ungrooved  static  ultimate 
tensile  strength  and  both  are  identical  for  ungrooved  specimens.  However, 
Juvinall' s  method  of  constructing  the  S-N  diagram  for  grooved  specimens 
differs  from  Shigley's  in  the  use  of  strength  reduction  modifiers. 

In  considering  the  development  of  the  S-N  diagram  for  grooved  specimens 
Juvinall  offers  a  more  detailed  account  of  notch-sensitivity  and  stress 
concentration  (7,  pp.  237-262] .  However,  his  modifications  for  estimating 
the  notched  endurance  strength  differs  only  slightly  from  Shigley's.  In 
the  absence  of  specific  data  Juvinall  recommends  that  the  grooved  endurance 
strength  be  estimated  from 


S 


e 


0.5  S  k 
u  e 


k 

a 


(3.3-11) 


where 


S  =  static  tensile  ultimate  strength 

k  =  load  constant 

e 

k,  =  size  factor  or  diameter  factor 

D 

k  =  surface  finish  factor 

a 

K£  =  fatigue  stress  concentration  factor. 

Juvinall  then  offers  the  following  modification  to  Eq.  (3.3-8),  the 
relationship  for  calculating  K£  from  q  and  K^: 

Kf  =  1  +  (Kt  -  1)  q  ka,  (3.3-12) 


-  ■frffltfriiiiinffrT-  iifa 


fcAii 
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and  explains  [7,  pp.  256-257]: 

"A  conservative  procedure  which  is  often  followed  in 
estimating  the  endurance  limit  of  notched  parts  with 
various  surface  finishes  is  to  (1)  multiply  the  basic 
endurance  limit  (0.5  Su)  by  the  surface  finish  factor 
k  and  (2)  divide  by  the  fatigue  stress-concentration 
factor  as  determined  from  Eq.  (3.3-8).  The 
surface  factor  should  of  course,  pertain  to  the  finish 
at  the  stress  raiser,  as  this  is  the  point  where  a 
fatigue  fracture  would  presumably  originate.  For  parts 
having  values  of  k  appreciably  less  than  1,  the  above 
procedure  probably  errs  somewhat  on  the  conservative  side. 

The  reason  is  that  surface  irregularities  reduce  the 
additional  damage  which  can  be  caused  by  a  notch  in  much 
the  same  way  as  internal  irregularities  do.  Equation 
(3.3-8)  reflects  the  influence  of  internal  irregularities 
in  reducing  the  severity  of  the  notch.  A  modified  version 
of  this  equation  is  suggested  which  uses  the  factor  k  to 
compensate  for  the  influence  of  surface  irregularities", 
in  the  form  of  Eq.  (3.3-12). 

In  discussing  the  load  factor,  kg,  [7,  pp.  226-231]  Juvinall 
suggests  that  in  reversed  bending  the  loading  really  has  no  effect  upon  the 
modified  endurance  strength,  thereby  recommending  that  be  equal  to  unity 
for  machine  parts  in  purely  reversed  bending.  Juvinall  then  suggests  the 
following  thumb  rules  for  determing  the  effect  of  loading  [7,  p.  231]: 

Reversed  or  rotating  bending: 

k  =  1.0. 
e 

Reversed  axial  loads: 

kg  =  0.9  with  no  bending. 

k  =  0.6  to  0.85  with  indeterminate  bending, 

e 
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Reversed  torsion: 

k  =  0.58  for  ductile  materials, 
e 

k  =0.80  for  cast  iron, 
e 


Juvinall  is  in  close  agreement  with  Shigley  in  estimating  the 
effect  of  the  size  factor,  k^.  He  provides  the  following  procedure  for 
estimating  [7,  p.  236]: 


k^  =  1.0  for  d  ^  0.4  in.,  and  for  either  bending  or 

torsion.  (3.3-15) 

k^  =  0.9  for  0.4  in.£  d  ^  2  in.  and  for  either  bending 

or  torsion.  (3.3-14) 

k^  =  1.0  for  axial  loading.  (3.3-15) 


For  estimating  the  stress  level  to  be  plotted  at  10°  cycles,  Juvinall 
introduces  a  modified  fatigue  stress-concentration  factor,  K^'  [7,  pp. 
260-262].  Juvinall  explains: 


"Thus  far  we  have  noticed  that  the  influence  of  stress 
raisers  on  strength  at  infinite  life  (or  a  very  large  number  of 
cycles)  is  represented  by  and  that  the  effect  of  stress 
raisers  on  static  strength  is  commonly  neglected  with  ductile 
materials.  Our  understanding  of  the  material  behavior  for  these 
two  limiting  cases  leads  us  to  expect  that  the  influence  of  a 
stress  raiser  would  increase  in  some  continuous  manner  over  the 
intermediate  or  finite-life  range.  This  is  indeed  the  case. 
Stress  concentration  factors  for  finite  life  are  designated  by 
K^,  where  1  £  Kf'  £  Kf." 


Juvinall  presents  a  curve  to  estimate  the  values  of  K^'  [7,  Fig.  13,26  , 

p.  260]  in  the  form  of  Fig. 3.3-5  ,  once  and  S  have  been  established. 

Juvinall  suggests  the  following  relationship  for  estimating  the  stress 

3 

level  at  10  cycles: 


0.9  S 


S103  cycles 


u 


K, 


(3.3-16) 
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Having  determined  the  necessary  parameters  for  using  Juvinalls* 
recommendations  the  S-N  diagram  may  now  be  plotted  for  the  same  specimen. 
Using  Juvinalls*  development  for  estimating  the  grooved  endurance  strength 
at  10°  cycles,  or  Eq.  (3.5-11),  we  have 


where 


and 


u 


0.5  S  k  k,  k  , 
u  e  b  a 


106,000  psi, 

1  [7,  p.  231] , 


kj,  =  1  [7,  p.  238], 


k  =  0.75  [7,  p.  234]. 

2 


Employing  q  =0.89  and  K  =  1.06,  Eq.  (3.3-12) 

v> 


gives 

Kf  =  1  +  (1.06  -  1) (0.89) (0.75), 
or 

K£  =  1.04. 

The  calculation  of  S  may  be  completed  by  substituting  K  into  Eq. (3.3-11) 

0  jt 

or 

„  (0.50) (106,000) (1)(1) (0.75) 

e  “  ’ 

e  1.04 

— i 

Sg  =  38.2  Kpsi. 

This  is  the  stress  level  at  the  10^  cycle  point  on  the  S-N  diagram. 
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Now  the  stress  level  at  the  10  cycle  point  must  be  calculated 
using  Eq.  (3.3-16).  Knowing  and  K^,  K^'  is  found  to  be  1.009  from 
Fig.  3.3-5.  Substituting  all  lvalues  into  Eq.  (3.3-16)  yields  the  10"5 
cycles  stress  level  as 


or 


S103 


cycles 


(0.9) (106,000) 
1.009 


^10^  cycles 


94.6  Kpsi. 


The  S-N  diagram  can  now  be  plotted,  as  shown  in  Fig.  3.3-6.  The 
corresponding  experimentally  determined  S-N  diagram  for  this  specimen 
is  superimposed  in  Fig.  3.3-6.  The  same  procedure  was  used  for  the 
same  specimen  but  for  a  groove  radius  of  r  =  0.125  in.  The  results  are 
superimposed  in  Fig.  3.3-6. 


COMPARISON  OF  JUVINALL'S  EMPIRICAL  METHOD  WITH  EXPERIMENTALLY 
DETERMINED  RESULTS 


Upon  comparing  in  Fig.  3.3-6  the  S-N  diagram  constructed  using 
Juvinall's  method  with  that  constructed  using  experimental  data 
exclusively,  it  is  found  that,  like  Shigley's  method,  Juvinall's  method 
gives  very  conservative  results.  This,  however,  is  not  surprising  con¬ 
sidering  the  close  proximity  of  the  two  methods.  The  mean  value  at  103 
cycles  by  Juvinall's  method  is  94.6  Kpsi,  whereas  the  mean  value  line 
on  the  experimentally  determined  S-N  diagram  crosses  the  ordinate  for 

3 

10  cycels  at  112  Kpsi.  Also  like  Shigley's  method,  Juvinall's  method 
gives  conservative  results  at  the  endurance  end  of  the  spectrum.  At 
106  cycles  Juvinall's  value  is  38.2  Kpsi  while  the  experimental 
endurance  value  is  54  Kpsi.  This  it  would  seem  that  Juvinall's  method 
becomes  more  conservative  toward  the  endurance  end  of  the  S-N  diagram. 
However,  when  there  is  a  lack  of  data,  or  time  to  run  experimental  tests, 
either  Juvinall's  or  Shigley's  method  can  be  used  and  all  errors  would 
be  on  the  "safe  side". 
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ALTERNATIVE  METHODS  FOR  CONSTRUCTING  S-N  DIAGRAMS  FOR  GROOVED  AND 
UNGROOVED  SPECIMENS 

Yokobori  [13,  pp.  195-336],  Madayag  [14,  pp.  110-113]  and  Osgood 
[15,  p.  371]  all  attempt  to  develop  equations  to  express  the  S-N  curvt-  in 
a  functional  form  for  both  ungrooved  and  grooved  specimens.  The  equations 
developed  by  all  three  authors  depend  upon  empirically  determined  constants. 
Yokobori  gives  the  following  equations  representing  attempts  to  express  the 


S-N  curve  in  a  functional  form  [13,  p.  195]: 

S  =  AN"y  ,  (3.3-17) 
S  =  a  -  b  log  N,  (3.3-18) 
log  (S  -  S  )  =  c  -  d  log  N,  (3.3-19) 
log  (S  -  Se)  =  e  -  f  log  (N  +  g).  (3.3-20) 


Equations  (3.3-17)  and  (3.3-18)apply  only  for  stresses  greater  than  the 
endurance  limit,  while  Eqs.  (3.3.19)and  (3.3-20)  include  the  endurance 
limit,  S  .  Due  to  the  wide  scatter  in  N  for  a  given  stress  level,  attempts 
to  determine  values  for  the  constants  in  Eqs.  (3.3-17)  through  (3.3-20) 
have  not  been  very  reliable  and  should  be  carried  out  on  a  statistical 
basis  [13,  p.  196].  It  should  be  pointed  out  that  Eqs.  (3.3-17)  through 
(3.3-20  )  can  be  used  for  both  grooved  and  ungrooved  parts  as  the  constants 
would  change  for  grooved  specimens. 

Madayag  [14,  pp.  110-113]  also  gives  equations  that  attempt  to  express 
the  S-N  curve  in  mathematical  form.  The  first  equation  is  that  proposed  by 
Weibull  [14,  p.113]  and  is  of  the  form 

(S  -  S. )(N  +  B)a  =  b,  (3.3-21) 

where  S  is  the  endurance  limit  and  a,  b  and  B  are  constants.  Valluri 
e 

[14,  p.  113]  proposed  the  following  equation  for  determining  N: 

2  log e(°u/°)  ioget(  a~  ai)/K3 
c  [(a  -  a.)/E]2  [(a  -  o’)/  a^2 


N 


(3.3-22) 
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where 


a  = 

maximum  cyclic  stress. 

o 1  = 

minimum  cyclic  stress. 

c? 

n 

ultimate  tensile  strength. 

ai= 

internal  stress. 

ii 

o 

material  constants. 

E  = 

Young's  modulus. 

Osgood  also  develops  a  relationship  for  the  S-N  curve  [15,  p.  371], 
For  ungTooved  steel  bars  in  rotating-bending,  Sg,  the  endurance  strength, 
was  derived  as, 

Sg  =  (1.0139)1/n(K)°*911,  (3.3-23) 

and  for  grooved  bars, 

Se  =  (1.0108)1/n(K)°’828,  (3.3-24) 

where  n  is  the  slope  of  the  true  stress-strain  curve  and  K  is  the  stress- 
intensity  factor. 

All  of  the  equations  given  above  depend  upon  empirically  determined 
constants.  Whether  or  not  realistic  values  for  these  constants  can  be 
determined,  it  is  obvious  that  the  methods  presented  by  Shigley  and 
Juvinall  are  much  more  convenient  for  the  designer  in  estimating  the  fatigue 
strength  of  a  given  material.  These  methods  only  depend  on  knowing  the 
tensile  strength  of  a  given  material  which  is  easily  obtained,  and  on 
applying  the  various  thumb  rules  and  modifying  factors,  to  come  up  with 
the  complete  S-N  diagram. 
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3.3.5  NOTCH  SENSITIVITY 

3. 3. 5.1  ASPECTS  OF  NOTCH  SENSITIVITY 

IN  FATIGUE 

In  discussing  fatigue  strength  and  thi  construction 
of  S-N  diagrams  in  Section  3.3.4  the  influence  of  the  various  stress 
raisers  such  as  surface  finish,  size  effects  and  stress  concentration 
was  presented.  Although  all  stress  raisers  should  be  accounted  for  in 
a  good  engineering  design  the  most  significant  and  most  difficult  to 
predict  is  that  due  to  stress  concentration. 

Stress  concentration  in  a  machine  part  may  be 
influenced  by  a  variety  of  factors  such  as  internal  cracks  or 
irregularities  caused  by  machining.  However,  the  most  common  cause 
of  stress  concentration  is  that  due  to  unavoidable  changes  in  the 

cross-sectional  area  of  the  machine  part,  such  as  a  hole,  groove,  or 
keyway.  In  order  to  account  for  the  decrease  in  strength  due  to  the 


different  materials  to  the  presence  of  a  groove  an  index,  q,  known  as 
notch  sensitivity  was  defined  in  Section  2. 2.-1  by  Eq.  (2.2-27), 
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q  = 


(2.2-27) 


Since  Kt  values  are  readily  available  the  notch  sensitivity  index  was 
defined  primarily  to  give  designers  a  method  to  calculate  K_  without 

T 

extensive  experimental  testing.  A  material  fully  notch  sensitive 
would  have  a  q  of  1  and  a  material  completely  insensitive  to  notches 
would  have  a  q  of  0.  However,  since  q  in  addition  to  material 
properties  is  also  affected  by  the  sharpness  of  the  groove,  the 
depth  of  the  groove,  etc.  very  little  reliable  notch  sensitivity 
data  has  been  generated  to  date  elsewhere.  The  problem  is  basically 
the  same  as  that  for  finding  Kfj  i.e.  extensive  testing  is  required 

to  determine  q.  Although  q  normally  assumes  values  between  zero  and 
1,  in  some  instances  values  of  q  greater  than  one  have  been  calculated. 
The  cause  of  this  may  be  traced  to  the  problem  of  maintaining  complete 


control  over  the  variability  involved  in  fatigue  testing.  Published 
data  by  NASA  researchers  illustrating  this  problem  is  depicted  in 
Fig.  3.3-7. 
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Table  3.3-6  -  Rotating  Beams  With  Circumferential  Grooves  (Cont'd). 
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Table  3.3-7  Check  Of  Computed  Strength-Reduction  Factor  (k1)  With 
Strength-Reduction  Factor  Determined  Directly  by 
Fatigue  Tests  (Gont ' d) . 
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Fig.  3.5-7  Notch  sensitivity  versus  notch  radius.  Vertical  lines 
represent  observed  "points".  Dashed  line  represents  ' 
Neuber's  theory  [16,  p.  47,  Fig.  24).  P 
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The  first  detailed  analysis  of  the  influence  of  stress 
gradients  and  material  micro-structure  on  the  values  of  stress  raisers 
was  made  by  Neuber  [19,  p.  17}.  Neuber  undertook  an  extensive  theo¬ 
retical  analysis  of  notch  effects,  and  developed  the  Neuber* s  technical 
stress-concentration  factor  designed  by  [19,  p.  16],  The  equation 
for  is  given  by 


K  =  1  + 
n 


(3.5-25) 


where  r  is  the  notch  radius  in  inches,  m  is  the  notch  flank  angle  in 
radians,  is  the  theoretical  geometric  stress  concentration  factor 
and  "a"  is  an  empirical  constant  which  represents  half  the  length  of 
an  equivalent  grain  of  the  material .  Values  of  "a"  are  found  by 
experiment,  to  be  portional  to  grain  size  in  a  very  general  way,  but 
are  actually  several  times  the  actual  size.  Kuhn  [17,  p.  27,  Fig.  5] 
has  published  curves  giving  values  for  various  materials,  which  is 
reproduced  in  Fig.  3.3-8.  This  curve,  however,  is  very  general  and 
there  still  is  much  discussion  and  research  involved  with  Neuber's 
material  constant. 

A  refinement  of  Neuber's  original  approach  was  developed  by 
Peterson  [5,  pp.  10-16].  This  refinement  is  convenient  to  use  and 
has  been  widely  accepted.  Peterson  bases  his  curves  on  the  notch 
sensitivity  factor,  q,  and  Eq.  3.3-8. 


■_  ]  -  |-  M  it  J^ux. * 
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Using  specimens  with  a  zero  flank  angle  he  reduces  Neuber's 
equation  to 


!  +  K  -  ! 


(3.3-26) 


Substitution  of  Eq.  (3.3-26)  into  Eq.  (3.3-8)  yields 


(3.3-27) 


Eq.  (3.3-27)  supplemented  with  values  of  "a"  from  Fig.  3.3-8  gives 
vaules  of  q  for  bending,  axial  and  torsional  loading.  Such  values 
are  plotted  in  Fig.  3.3-9.  Most  of  the  notch  sensitivity  work  under¬ 
taken  in  this  research  is  based  on  Peterson's  refinement  of  Neuber's 
theories.  Figs.  3.3-10  and  3.3-11,  also  generated  by  Peterson, 
relate  q  to  geometric  configuration  and  are  based  on  Neubex's  work. 

The  curves  are  completely  theoretical  and  do  not  directly  involve 
Neuber's  material  constant.  The  set  of  equations  derived  by  Neuber 
for  K  and  used  by  Peterson  to  plot  his  K  curves  given  in  Figs.  2.2-7, 
2.2-8,  2.2-11,  and  2.2-12  for  circumferentially  grooved,  round  rods 
in  bending  are  given  in  Fig.  3.3-12. 

It  has  been  shown  by  many  test  that  use  of  the  full  theoretical 
value  of  Kt  in  fatigue  design  will  generally  give  results  on  the 
conservative  side,  but  for  materials  of  any  reasonable  degree  of 
ductility  plastic  strain  occurs  in  the  notches  and  reduces  the  stress 
concentration  effect  to  a  value  some./hat  below  that  calculated  from 
pure  geometry.  This,  leads  into  another  definition  of  notch  sensi¬ 
tivity  in  terms  of  and  K^.  It  has  been  found  experimentally  that 
there  is  a  limiting  small  value  of  the  notch  radius  below  which  there 
is  no  additional  stress  concentration  effect.  The  ratio  between  the 
apparent  increase  in  local  stress  in  fatigue  and  the  increase  predicted 
by  the  elastic  theory  of  stress  concentration  is  another  way  of 
explaining  the  parameter,  q  [15,  pp.  103-104]. 


Fig.  3.3-9  Notch  sensitivity  values  versus  notch  radius  [7,  p.  255, 
Fig.  13.21]. 


3. 3. 5. 2  COMPARISON  OF  THEORETICAL  AND  ONR  EXPERIMENTAL 


154 


NOTCH  SENSITIVITY  RESULTS 


Notch  sensitivity  is  a  highly  sensitive  parameter  and  varies 
significantly  depending  upon  (1)  the  material,  (2)  the  type  and 
severity  of  the  notch,  and  (3)  the  type  and  severity  of  loading 
[15.  p-  104  J  .  For  this  research  fully  reversed  bending  is  the  only 
loading  considered.  It  has  been  observed,  as  in  Table  3.3-8,  that 
some  alloys  show  a  equal  to  one  indicating  nc  reduction  in  fatigue 

strength  due  a  notch  while  in  this  research  some  configurations  of 
AISI  1038  steel  have  a  K,.  value  as  high  as  1.54,  as  given  in  Table 

3.3- 9,  indicating  a  high  reduction  in  strength  due  to  the  presence 
of  the  notch. 

The  results  of  this  research  for  q  and  Kf  are  given  in  Table 

3.3- 9  and  are  plotted  in  Figs.  3.3-13  and  3.3-14  and  confirm  the 
accepted  theory  as  illustrated  in  Figs.  3.3-15,  3.3-16,  and  3.3-17. 

However,  it  should  be  noted  that  and  q  are  highly 
sensitive  parameters  and  great  care  must  be  taken  in  specimen 
preparation  and  testing. 
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Fatigue  Stress  Concentration  Factor,  Kf,  vs.  r/D  for  a  grooved  shaft  in  bending  AISI  1038  steel  notch 
rod,  endurance  at  3  x  106  cycles,  rod  diameter  D  =  0.375  in.,  base  diameter,  d  =  0.2700  in.  varied/radii . 


NOTCH  RADIUS,  f,  IN. 


3-15  Notch  sensitivity  versus  notch  radius 
[18,  p.  526,  Fig.  4] 
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The  AISI  1018,  AISI  4340  and  AISI  4130  are  under 
investigation  now.  Their  notch  sensivity  and  fatigue  notch 
stress -concentration  factors  will  be  presented  in  the  next 
report. 
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3. 3. 5. 3  RESULTS  DRAWN  PROM  NOTCH  SENSITIVITY  TESTING 

The  results  of  these  notch  sensitivity  tests,  summarized  in 
Table  3.3-9,  are  based  upon  data  presented  in  Section  II  of  this 
report.  The  objective  was  to  determine  the  effect  of  decreasing 
notch  radius  upon  fatigue  stress  concentration  factor,  Kf,  and  notch 

sensitivity,  q. 

The  response  to  notch  configuration  plotted  in  Figs.  3.3-13 
and  3.3-14  indicate:*  that  a  decrease  in  Kf  from  1.540  to  1.155  and  an 
increase  in  q  from  0.636  to  0.968  results  from  an  increase  in  the  r/d 
ratio.  Although  the  fatigue  stress -concentration  factor,  Kf,  increases 

as  the  r/d  ratio  decreases  the  notch  sensitivity  factor,  q,  decreases. 
This  indicates  that  as  r/d  decreases  the  material  becomes  less  notch 
sensitive  although  the  strength  of  the  specimen  is  decreased  due  to 
the  presence  of  the  notch. 

The  results  presented  in  Table  3.3-9  and  the  curves  plotted 
in  Figs.  3.3-13  and  3.3-14  may  be  employed  for  the  design  of  shafts 
of  similar  configurations  in  reversed  bending. 

3.4  AXIAL  FATIGUE  RESEARCH  MACHINE 

3.4.1  TENSILE  DISTRIBUTIONAL  PROPERTIES 

Table  8.3-1  presents  the  experimental  data  generated  to 
determine  the  distributional  static  ultimate,  yield,  and  breaking 
strengths  of  AISI1018  steel  rod  used  in  this  research.  Table  8.3-2 
presents  the  hardness  and  percent  elongation  for  the  static  test 
specimens. 
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3.4.2  DISTRIBUTIONAL  GOODMAN  DIAGRAMS  FOR  AISI 
1018  STEEL 

3. 4. 2.1  DETERMINATION  OF  THE  STRENGTH  PARAMETERS 
AT  SPECIFIC  STRESS  RATIOS 

Fatigue  tests  were  conducted  using  AISI  1018  steel 
specimens  -to  determine  the  endurance  strength  at  2.0  x  10&  cycles  of 
life.  The  specimens  were  subjected  to  alternating  axial  stresses 
superimposed  onto  mean  tension  stresses  at  alternating-to-mean  stress 
ratios,  rs,  of  °°,  2.0,  1.0,  0.4  and  0.2. 

The  staircase  test  procedure  discussed  in  Section 

3.2.3  was  followed  and  the  test  results  were  reduced  using  the 
standardized  data  reduction  procedure  presented  in  that  section. 

A  PDP-8  computer  program  was  developed  and  used  to  calculate  the  mean 
and  standard  deviation  of  the  strength-to-failure  (endurance)  data 
and  is  presented  in  Appendix  E. 

Once  the  mean  and  standard  deviation  of  the  stress 
vector,  Sy,  is  calculated;  the  mean  and  ‘  3  u  envelopes  for  each 

stress  ratio  can  be  located  and  a  distributional  Goodman  diagram  for 
specimens  of  a  specific  material  and  geometry  may  be  constructed. 

3. 4. 2. 2  CONVERSION  OF  STAIRCASE  TEST  RESULTS  INTO 

DISTRIBUTIONAL  GOODMAN  DIAGRAMS 

Goodman  diagrams  are  plotted  with  the  mean  stress 
located  on  the  abscissa  and  the  alternating  stress  on  the  ordinate. 
The  results  of  the  staircase  test  provides  the  mean  and  standard 


164 


deviation  of  the  stress  vector,  Sv.  Thereby  s Hewing  the  stress 

vector  distribution  to  be  plotted  directly;  with  the  stress  ratio, 
rs,  determining  the  angle,  a,  that  the  stress  vector  makes  in 

relation  to  the  abscissa,  as  shown  in  Fig.  3.4-1. 

The  curve  drawn  through  the  mean  value,  Sy,  of  each 

stress  vector  distribution  at  the  various  stress  ratios  constitutes 
the  mean  of  the  Goodman  diagram.  Similarly,  the  curves  connecting 
the  Sv  +  3  a  and  Sv  -  3  a  values  provide  the  ±  3  cr  envelopes,  hence 

the  distributional  Goodman  diagram. 

3. 4. 2. 3  COMPARISON  OF  EXPERIMENTAL  FATIGUE  TEST  RESULTS 
WITH  ESTABLISHED  FAILURE  THEORIES 

The  Goodman  diagram  constructed  from  experimental 
results  is  compared  with  curves  drawn  theoretically  using  the 
Distortion  Energy  Theory  and  the  Maximum  Shear  Stress  Theory.  The 
Goodman  diagram  constructed  in  accordance  with  the  Distortion  Energy 
Theory  is  described  by  the  von  Mises-Hencky  ellipse  and  is  defined 
by 


The  Goodman  diagram  constructed  in  accordance  to  the  Maximum  Shear 
Stress  Theory  is  slightly  more  conservative  than  the  von  Mises-Hencky 
ellipse  and  is  known  as  the  Gerber  parabola  defined  by 


le  of  plot  of  a  stress  vector  distribution  for 
ess  ratio,  rs  =  1. 
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Sa 

Se 


2 


=  1  . 


(3.4-2) 


The  curves  generated  by  using  the  established 
failure  theories  are  plotted  by  calculating  the  stress  vector,  Sv, 

corresponding  to  each  stress  ratio  and  using  the  established 
endurance  strength,  Se,  and  the  static  ultimate  tensile  strength,  Su, 
for  the  specified  material.  For  a  selected  mean  stress,  Sm,  and 
stress  ratio,  rs,  the  stress  vector  can  be  computed  from 

Sv  =  Sm  (1  +  xs)h  ,  (3.4-3) 


where  Sm  is  calculated  using  Eq.  (3.4-1)  or  Eq.  (3.4-2)  depending 

upon  the  failure  theory  selected. 

Points  on  the  von  Mises-Hencky  ellipse  are 
computed  by  selecting  several  stress  ratios  and  substituting  the 
appropriate  Se  and  Su  values  into  the  following  equation: 

S  -  S 

u  e  ^ 

Sme  (Se2  +  rs2  Su2)^ 

The  corresponding  stress  vector,  Sv,  value  is  calculated  using 
Eq.  (3.4-3)  and  substituting  the  value  of  Sm  found  from  Eq.  (3.4-4) 
into  Eq.  (3.4-3). 

If  the  Gerber  parabola  is  used  to  describe  the 
Goodman  diagram  the  following  relationship  is  used  to  calculate 


( 
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S  = 
m 
P 


2  4  2  2  2 

(rs  Su  +  4  Su  •  Se  )  -  rs  •  Su 

2  S. 


(3.4-S) 


Equation  (3.4-3)  is  then  used  to  calculate  the  stress  vector,  Sy, 
corresponding  to  the  Sra  value  found  by  Eq.  (3.4-5). 

A  PDP-8  computer  program  has  been  developed  which 
identifies  the  points  on  the  von  Mises-Hencky  ellipse  and  the  Gerber 
parabola  when  Se  and  Su  are  known  and  values  of  rs  are  specified. 

The  program  is  presented  in  Appendix  F  . 

The  curves  generated  by  the  von  Mises-Hencky 
ellipse,  the  Gerber  parabola  and  the  ONR  experimental  test  results 
for  AISI  1018  steel  are  shown  in  Fig.  (3.4-2).  It  is  apparent  that 
neither  the  ellipse  nor  the  parabola  fit  the  experimental  data, 
because  the  experimental  data  curve  lies  between  these  two. 


3.4. 2.4  METHODOLOGY  FOR  FINDING  THE  BEST-FIT  CURVE 


THROUGH  THE  EXPERIMENTAL  DATA  POINTS 


A  comparison  of  Eq.  (3,4-1),  describing  the 
von  Mises-Hencky  ellipse,  and  Eq.  (3.4-2),  describing  the  Gerber 
parabola,  reveals  that  both  theories  can  be  expressed  generally 
in  the  same  form  by 


sa> 


•(*)’■■■ 


(3.4-6) 


where  a  =  1  for  the  parabola  and  a  =  2  for  the  ellipse.  Since  the 


T 

i 


experimental  data  for  AISI  1018  steel  lie  between  the  ellipse  and  the 
parabola,  there  exists  a  value  of  the  exponent,  a,  which  will  describe 
the  experimental  curve.  Transformation' of  Eq.  (3.4-6)  into  a  linear 
relationship  allows  the  use  of  the  standard  "least  squares"  method  of 
curve  fitting  to  solve  for  the  value  of  the  exponent,  a. 

Substitution  of  y  =  — — •  and  x  =  Sm/Su  into 

be 

Eq.  (3.4-6)  yields 


a 

y 


Taking  the  natural  logarithm  of  both  sides  of  the 
transformed  Eq.  (3.4-6)  yields 

a  In  y  =  in  (1  -  x2)  . 

Further  substitution  of  Y  =  In  y  and  X  =  In  (1  -  *2,  int„  Eq.  (3.4-7j 
gives 


a  Y  =  X  , 


or 


Y  = 


(3.4-8) 


I 

i 
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Equation  3.4-8  is  now  the  linear  form  of  Eq.  (3.4-6)  and  is  in  the 


form  of  a  straight  line  where  the  y  -  intercept  is  zero  and  the  slope 
is  1/a. 

A  CDC-6400  FORTRAN  computer  program  has  been 
developed  which  calculates  the  slope,  l/a>  from  experimental  data. 
Inputs  to  the  program  include  the  values  of  Sa  and  Su  for  the 

material  and  values  of  Sv  and  rs  from  the  staircase  testing  program. 

Outputs  of  the  program  include  the  exponent,  a,  for  Eq.  (3.4-6)  and 
the  correlation  coefficient  for  the  fit  of  the  line  describing  the 
data  points.  Additional  outputs  of  the  program  include  any  specified 
number  of  Sv  and  rs  coordinates  to  provide  for  the  plotting  of  the 
experimental  curve. 

Analysis  of  the  program  results  to  date  have 
revealed  two  significant  points: 

1.  The  exponent  calculated  by  the  program  is 
sensitive  to  the  presence  or  absence  of 
data  in  the  stress  ratio  range  of  0.2  to 
0.4. 

2.  Inclusion  of  the  end  point  where  Ys  =  0 
resulting  in  Sm  =  Su  cannot  be  included 
directly  because  it  results  in  taking  the 
natural  logarithim  of  zero  which  is 
undefined.  However,  inclusion  of  the  end 
point  is  needed  to  assure  the  proper 
curvature  of  the  data  line.  This  problem 
was  circumvented  by  making  use  of  the 


extreme  accuracy  of  the  CDC-6400  computer  and 
offsetting  the  end  point,  by  0.0001  psi  and 

the  stress  ratio,  rs,  by  0.0001. 

The  significance  of  the  first  point  lies  in  the 
fact  that  the  difference  between  the  ellipse  and  the  parabola  is  the 
maximum  in  the  same  range  of  stress  ratios  of  0.2  to  0.4. 

Consequently  more  fatigue  tests  must  be  conducted  in  this  range  of 
stress  ratios  to  determine  the  correct  exponent  of  Eq.  (3.4-6). 

Figure  3.4-2  is,  therefore,  tentative.  Fatigue  tests  for  AISI  1018 
steel  are  continuing  in  this  stress  ratio  range.  The  results  will 
be  included  in  the  next  report. 

The  CDC-6400  program  briefly  described  above  is 
being  expanded  to  include  the  calculation  of  data  points  for  the 
dis tributional  Goodman  diagrams  i  3  a  envelopes.  The  final  program 
will  be  included  in  the  next  report. 
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4.0  OVERALL  CONCLUSIONS 

4.1  WIRE  FATIGUE  RESEARCH  MACHINE  TESTS 

1.  Static  tensile  tests  were  successfully  completed  on  0.040 
in.  diameter  wire  specimens  made  of  AISI  1018,  1038,  4130 
and  4340  steel.  The  distributional  strength  parameters 
are  given  in  Table  3.2-1. 

2.  The  ultimate  strength  mean  values  for  wires  are  in  general 
agreement  with  conventional  tabled  values  for  the  specified 
material. 

3.  The  coefficients  of  variation  vary  inversely  with  the  values 
of  ultimate  strength.  The  coefficient  of  variation  values 
range  from  11%  for  AISI  4340  steel  to  25%  for  AISI  1038 
steel,  values  somewhat  higher  than  for  steel  rods  of  the 
same  material. 

4.  Static  strength  distributions  can  be  effectively  represented 
by  the  normal  distribution. 

5.  Cycles-to-failure  tests  were  successfully  completed  on  the 
wires  of  the  same  materials  and  diameters.  Distributional 
parameters  calculated  for  the  normal,  log  normal,  and 
Weibull  distributions  are  listed  in  Tables  3.2-2  through 
3.2-5. 

6.  Preparation  of  distributional  S-N  diagrams  have  not  been 
completed  because  endurance  testing  is  still  in  progress. 
Preliminary  plots  shown  in  Figs.  3.2-1  through  3.2-4 
confirm  that  distributional  means  and  envelopes  can  be 


Preceding  page  blank 
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approximated  by  straight  lines  when  the  logarithms  of  cycles 
to  failure  are  plotted  vs  the  logarithm  of  stress. 

7.  The  cycles-to-failure  distributions  have  large  variances 
with  coefficients  of  variation  values  as  high  as  50%  at  the 
lower  stress  levels.  This  variability  further  emphasizes 
the  importance  of  design  by  reliability  techniques  to 
replace  conventional  design  based  on  estimates  of  determi¬ 
nistic  values. 

8.  Based  on  the  results  of  the  K-S  and  Chi-Squared  goodness- 
of-fit  tests,  the  coefficients  of  skewness  and  kurtosis, 
and  phenomenological  considerations  the  lognormal  distri¬ 
bution  provides  the  best  fit  to  the  cycles-to-failure  data. 
However,  the  normal  distribution  also  appears  to  provide  a 
useful  fit. 

9.  The  three-parameter  Weibull  distribution  can  be  considered 
to  provide  an  acceptable  fit  for  cycles-to-failure  data 
based  on  the  K-S  goodness-of-fit  test  results.  As  the 
Chi-Squared  test  is  more  discriminating  for  small  sample 
sizes,  its  result  show  that  the  Weibull  distribution  should 
be  used  with  care  for  sample  sizes  of  35  or  fewer. 
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WIEDEMANN  RESEARCH  MACHINE  TESTS 

1.  Static  tensile  test  results  cn  AISI  1038,  4130  and  4340 

steel  ungrooved  rods  of  0.2S-in.  nominal  diameter  are 

summarized  in  Table  3.3-1,  The  following  conclusions  may 

b3  drawn: 

a)  Ultimate  strength  mean  values  are  in  general  agreement 
with  published  properties. 

b)  The  coefficient  of  variation  for  the  ultimate  strength 
varies  from  1%  to  2%  which  attests  to  the  homogeneity 
of  the  materials  and  the  -high  precision  in  the  tests. 

c)  The  percent  elongation  of  the  AISI  1038  steel  specimens 
had  a  mean  of  20.9%  as  compared  with  the  mean  of  8.5% 
for  AISI  4130  steel  specimens.  These  results  are  as 
expected  since  AISI  1038  steel  is  considerably  softer 
than  AISI  4130  steel. 

d)  Because  elongation  is  accompanied  by  a  decrease  in  the 
cross-sectional  area,  the  mean  breaking  strength  of  the 
AISI  1038  steel  specimens,  based  on  final  diameters,  is 
larger  than  for  the  AISI  4130  specimens.  These  results 
illustrate  that  caution  must  be  excercised  when  using 
the  breaking  strength  as  a  design  criterion. 

2.  Cycles-to-failure  test  results  for  AISI  4130  steel  rod 
specimens  with  0,0250"  and  0.125*'  groove  radii  are  complied 
in  Tables  3.3-3  and  3.3-3,  respecitvely.  The  following 
conclusions  may  be  drawn: 
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a')  S-N  diagrams  plotted  on  log-log  scales  using  normal 
cycles-to-failure  distribution  parameters  confirmed 
that  distributional  S-N  diagrams  approximated  by 
strai ght-line  loci  for  the  mean  and  ±  3  a  envelopes  are 
good  representations  of  the  data. 

b)  The  variability  of  the  cycles  to  failure  increased  as 
the  stress  levels  decreased  so  that  the  distributional 
S-N  diagram  retained  the  straight  line  fit  with  widening 
±  3 a  envelopes  as  the  stress  levels  decreased. 

c)  A  comparison  of  the  experimentally  generated  S-N 
diagrams  witr>  existing  empirical  methodologies  showed 
that  (1)  the  empirical  methodologies  are  conservative 
and  (2)  the  distributional  information  provided  by  the 
experimental  diagrams  allows  for  improvements  in  design 
accuracy. 

3.  Endurance  tests  on  AISI  1038  steel  rod  specimens  with 
various  groove  radii  under  reversed  bending  stresses  were 
conducted  using  the  staircase  methodology.  The  following 
conclusions  may  be  drawn  from  the  results  tabulated  in 
Table  3.3-4, 

a)  Specimens  having  a  groove  radius  of  1.87"  which  provides 
a  stress  concentration  factor  of  approximately  1.0, 
have  a  mean  endurance  strength  of  37,100  psi.  As  the 
groove  radius  is  decreased  to  0.031",  the  mean  endurance 
strength  decreases  to  23,800  psi. 
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b)  A  plot  of  the  endurance  strength  for  different  groove 
radii  provided  for  an  analysis  of  notch  sensitivity. 

The  values  of  notch  sensitivity  given  in  Table  3.3-9 
agreed  very  closely  with  published  results. 

c)  When  conducting  staircase  endurance  tests,  it  was  found 

that  to  minimize  invalid  test  points  the  initial  stress 

level  should  be  0.54  times  the  ultimate  strength  (S^) 

for  groove  radii  of  1.87  or  larger.  For  smaller  groove 

radii  the  initial  stress  level  should  be  0.54  S  /K... 

u  f 

4.3  AXIAL  FATIGUE  RESEARCH  MACHINE  TESTS 

Endurance  staircase  testing  was  completed  on  AISI  1018  steel 
specimens  at  alternating  to  mean  stress  ratios  of  2.0,  and  1.0, 
and  is  in  progress  for  stress  ratios  of  0.4  and  0.2.  The  following 
conclusions  may  be  drawn  from  the  results  given  in  Table  2.3-1  and 
Figs.  9.3-1  thru  9.3-3. 

1.  At  the  stress  ratio  of  »,  the  endurance  strength  was 
0.434  Su  which  is  in  general  agreement  with  prior  results 
obtained  for  AISI  1038,  4130,  and  4340  steels.  Thus 

0.45  Su  provides  a  good  starting  stress  level  for  staircase 
type  endurance  tests  with  mean  and  alternating  axial  loads. 

2.  While  the  experimental  Goodman  diagram  more  closely 
approximates  the  Gerber  parabola  than  the  ellipse,  the 
correct  curve  is  somewhere  bewteen  the  two.  A  least 
square  regression  analysis  is  being  developed  to  determine 
the  best  representation  of  the  Goodman  diagram  data. 


3.  The  largest  difference  between  the  Gerber  parabola  and  the 
von  Mises-Hencky  ellipse  occurs  in  the  stress  ratio  range 


of  0.2  to  0.4;  consequently  more  experimental  data  must 
be  obtained  within  this  range  of  stress  ratios  to  assure 
a  more  accurate  determination  of  the  Goodman  diagrams  for 
their  more  effective  use  in  design. 
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5.0  RECOMMENDATIONS 

1*  Experimental  research  in  progress  on  AISI  1018,  1038,  4130,  and 
4340  steel  specimens  should  be  completed.  Priority  should  be 
given  to  the  following  incomplete  efforts: 

a.  Development  of  distributional  S-N  diagrams  for  0.040  diameter 
wire  specimens  of  the  above  four  steels. 

b.  Development  Oi  the  distributional  Goodman  diagram  for  AISI  2018 
steel  rods  subjected  to  axial  alternating  and  mean  loads. 
Experimental  fatigue  data  is  needed  for  stress  ratios  of  0.4 

and  0.2  to  assure  accurate  determination  of  the  curvature  of 
the  Goodman  surface. 

c.  Development  of  distributional  S-N  diagrams  for  AISI  4150  0.0937" 
diameter  steel  rods  with  groove  radii  of  0.031"  and  0.062", 
subjected  to  reversed  bending. 

d.  Completion  of  endurance  strength  staircase  testing  of  AISI  1038 
and  4150  steel  rods  for  notch  sensitivity  analysis. 

2.  Experimental  data  generated  under  this  research  effort  has  been 
reported  in  a  series  of  technical  reports.  Distributional 
data  tables  should  be  developed  and  published  which  present 
experimentally  validated  data  for  use  by  designers  and  reliability 
engineers.  The  data  tables  should  also  include  published  data 
available  from  other  sources. 


Analysis  of  the  data  tables  should  be  performed  and  generalized 
thumb  rules  for  the  use  of  the  data  should  be  developed.  The 
thumb  rules  should  not  only  provide  assistance  in  the  use  of  the 
data  tables  but  should  include  rules  for  extending  the  data  tc 
other  materials  through  appropriate  relationships. 

A  combination  of  experimental  research  and  theoretical  study 
should  be  completed  to  determine  the  most  accurate  failure 
governing  criteria  for  AISI  1018,  1038,  4130  and  4340  steels 
subjected  to  fatigue  loading.  The  following  are  recommended: 

a)  Determine  functional  relationships  among  the  fatigue 
failure  governing  criteria,  type  of  load,  notch  sensitivity, 
and  size. 

b)  Analyze  singularly  and  jointly  the  distributions  of  notch 
sensitivity  and  .size  effects,  and  synthesize  these  distri¬ 
butions  into  failure  governing  stress  and  strength 
distributions . 

c)  Determine  the  functional  relationships  by  regression 
analysis  between  fatigue  strength  and  the  notch  sensitivity 
factor,  and  fatigue  strength  and  the  size  factor. 

The  effect  on  service  life  and  reliability  of  rotating  components 
of  cumulative  fatigue  loads  should  be  determined. 
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